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SUMMARY 


It has been known since the late tw;enties (Piercy and Richard- 
son) that unexpectedly high velocity and presumably pressure fluctua- 
tions are found in the stagnation region of blunt-nosed bodies placed in 
wind tunnels with uniform streams. Very recently, visualization techniques 
have revealed frequent formations of streamwise vortices both steady and 
vinsteady, These fluctuations and vortices are known to increase sub- 
stantially the heat and mass transfer in the stagnation region when 
r.m.s, free-stream turbulence exceeds approximately 1% of free-stream 
speed, U , and may also cause increased acoustic radiation. The 
boundary layer, though perturbed by the streamwise vortices tends to 
remain laminar in the region of accelerations around the head of the 
bluff bodies. However, these vortical disturbances are likely to con- 
tribute to earlier-than-expected transition to txirbulence farther 
downstream so that prevention of such "seeding" should delay transition. 

In an attempt to understand the phenomena behind these impor- 
tant practical flow manifestations, a critical review of the theoreti- 
cal and experimental evidence was made . Current theory is revealed to 
be incomplete, incorrect, or inapplicable to the phenomena observed 
experimentally. The formalistic approach via the "principle" of 
exchange of instabilities should most likely be replaced by a "forced 
disturbance" approach. Also, many false conclusions were reached by 
ignoring the fact that treatment of the base and perturbed flows in Hiemenz 
coordinate >9 is asymptotic in nature. Almost surely the techniques 
of matched asymptotic expansions will have to be used to capture 
correctly the diffusive and vorticity amplifying processes of the 
disturbances vis-a-vis the mean-flow boundary layer and outer poten- 
tial field as tf and y/diameter approach infinity. There are also 
serious uncertainties in the experiments which are discussed in 
detail. 
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1 INTRODUCTION 


The purpose of this essay is to summarize and illustrate ej^jerimental 
informa cion on flows in proximity of stagnation regions of cylindrical 
bodies and to examine critically the notions of flow instability, primarily 
as they might apply to these flows. The a pproacli Is expository , with 
emphasis on the structure of the phenomena and of the associated equations. 

The approach is also a searching one, because as new facts are iincovered, 
new questions also arise, especially with respect to our idealization of 
the phenomena . 

One might call the stagnation flows singular for they are physically 
complex, e.g. exhibiting singularly increased heat transfer response to 
even moderate free-stream turbulence, and also because their analyses in 
effect rely on singular perturbations. A glance at Figs .5, 7, 8, 11 and 12,’ 
samples of the most organized features of these flows, will convey an 
introductory impression of their complexity. 

The immediate motivation for focusing on these flows was to understand 
more about mechanisms which could bring about commonly observed streamwise 
vorticity in nominally two-dimensional boundary layers, Morkovin (1979). 

However, this report will not dwell on the facts and scientific or 
technological consequences of the occurrence of such strearawj.se vortical 
flows downstream of blunt (or sharp) leading edges nor on the rather 
controversial information and theories of the enhanced heat transfer. It 
will focus primarily on the curious physical and mathematical aspects of these 

unsteady, vortical, three-dimensional flow fields in stagnation regions. The 
rather discordant information on the enhanced heat transfer, found for instance 
in Kestin (1966) , Smith and Kuethe (1966) , Kayalar (1969) , Lowery and Vachon (1975) , 
(1975) , Miyazaki and Sparrow (1977) , Hoshizaki et al (1975) , Swigart (1977) 
and others , will be brought in only insofar as it provides evidence for 
specific flow behavior. 

There is also much disagreement on the character of stagnation flow 
fields themselves and especially on their theoretical interpretations and 
modeling. In order to gain additional perspective, we shall need to 
consider basic questions of instability phenomena and their mathematical 
formulation. Relevant issues of temporal versus spatial stability, of 
receptivity to free-stream disturbances, of its relation to initial value 
problems of instability and of nonlinearity are discussed in Appendix 1. 



However, only the needed, still tentative conclusions are utilized in the 
body of the report. 

The initial free-stream disturbances convected toward the cylinders 
are generally turbulent but the preseparation boundary layers on the 
cylinders, while displaying unsteady and random features, do not exhibit 
basic characteristics of turbulent wall layers (even when enhanced heat 
transfer is observed) . Here, and for stability problems in general, we 
need to sharpen our concepts and terminology of what constitutes turbulent 
flows. Appendix 3 proposes operational definitions of turbulence and 
provides some background for considerations of a class of quasi- laminar 
flows , such as the boundary layers under study here , which could be called 
buffeted laminar wall layers . 

Some of the experimental and theoretical information discussed in this 
essay is at least partially contradictory, not an uncommon situation in 
problems dealing with instabilities, transition to turbulence and 
turbulence itself. The author presents here points of view close to a 
perceived consensus of specialists with whom he had discussed the issues 
over the last three years. Contradictory and unconfirmed results have 
been recognized as a major block in the progress of transition research so 
that the U.S. Boundary Layer Transition Study Group adopted a set of 
guidelines to remedy such uncertainties and avoid miscues for further 
research , Reshotko (1976), p345. Guideline No. 4 states: "Whenever possible, 

tests should involve more than one facility; tests should have ranges of 
overlapping parameters, and whenever possible, experiments should have 
redundancy in transition measurements". It has been broadened to theoretical 
research where by a "facility" we understand a theoretical model with its 
computer program", Morkovin (1978). We shall have to invoke this generalized 
attitude , designated henceforth as Gen. Reshotko 4 , time and again in order 
to convey proper qualifications on the degree of reliability of various cited 
results and judgements. Despite this caution, it can happen that future 
experiments and theories will modify significantly the currently perceived 
consensus as here interpreted. In fact, a major purpose of describing 
the new perspective on the stagnation flows is to stimulate further theory 
and experiments specifically designed to test different elements of this 
perspective in accordance with Gen. Reshotko 4. 
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The exposition in Section 2 starts with an account of the piecemeal 
information on the steady or mean viscous flows around cylindrical objects. 

In particular the stagnation Hiemenz flow, lacking any characteristic lengths, 
and its imbedding within outer irrotational flows are described. It is anti- 
cipated that an even more severe problem of rigorous matching of inner and 
outer flows must be faced if the manner in which free-stream inhomogeneities 
and disturbances bring about imstable responses in the stagnation region is 
to be understood. In Section 3 the different characteristics and specifica- 
tions of these disturbances are first examined. An account of experimental 
observations of steady vortex formations upstream of cylindrical obstacles 
follows with an attempt at interpretation in terms of prototypes of nonlinear 
stability theory. The vortical instabilities in the stagnation region are 
commonly linked with the Taylor vortex instability in annular Couette flow and 

II , , 

with the Gortler streamwise-vortex xnstability in concave boundary layers , In 
order to gain further insight into instability phenomena in general, the 
differences in the physical and mathematical structure of these three super- 
ficially similar systems are outlined and discussed in Section 4. Considera- 
tion of disturbances and forced responses of the systems appears necessary for 
deeper understanding; see also Appendix 1. 

In Section 5, various processes and mechanisms are considered and 
illustrated through linearized equations , The Gortler-Hanmerlin perturbation 
model is examined with emphasis on its possible limitations in temporal depen- 
dence and in its confinement to the inner domain which lacks a characteristic 
length. Section 6 contains a careful account of the more detailed stagnation- 
flow experiments with vmsteadiness and turbulence in free stream. An overview 
of the important issues and the outlook for further progress is contained in 
Section 7. Useful information on various aspects of Gortler 's linear model is 
given in Appendix 4. 

The rather unconventional views presented here in a searching spirit 
were first written up in May 1976, Eighteen copies of that report were circu- 
lated in the United States, Europe, and Japan to elicit critique and sugges- 
tions and to provide an opportunity for private prepublication rebuttals for 
those whose earlier work was questioned in the text. The present report 
incorporates, then, not only the information and views from publications since 
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1976 but also the advice and comments of a number of discussors^ in particular' 
Profs. G, Inger, E, Roshotko, and I, Tani. Figures based on the reviewed 
piablications are explicitly credited, The original NASA study under Grant 
NSG-1120 was guided^ with understanding, by Ivan Beckwith, The 1979 revision 
was supported under AFOSR Contract F49620^77-c;'^013 concerned with instabili- 
ties and transition to turbulence. 

Use of trade names or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration, 
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2 CHARACTERISTICS OF THE OVERALL FLOW 
a. Piecemeal Modeling of Viscous Flows around Bodies, Matching eind 
Patching 

Before we can study the behavior of disturbances in a flow field we 
expect to have adequate knowledge of the undisturbed base flow . For 
slightly viscous flows around bluff cylindrical bodies the theoretical 
and experimental knowledge is piecemeal, resting on a series of idea- 
lizations and characterizations of local behavior. It is instructive 
for the rationale of handling the perturbed flow to examine the usual 
approach to the solution of the underlying Navier-Stokes equations for 
the base flow at Reynolds numbers in excess of say 1000, The Navier- 
Stokes equations are solved as such only in the immediate vicinity of 
the front stagnation line, where the classical solution of Hiemenz 
(Schlichting 1968, pp 87-90) demonstrates the existence of a thin viscous 
boundary layer without assuming its existence a priori. To continue the 
Hiemenz solution beyond the stagnation region, approximate boundary layer 
equations with an assumed outer potential flow are solved up to an ill- 
defined separation line where the wake begins (e.g. Schlichting 1968, pp 
154-162) . There is no satisfactory Navier-Stokes solution for the wake 
and rather crude semiempirical modeling is resorted to for assessments 
of its significant steady and unsteady influence on the potential flow 
around and upstream of the cylinder, (Only Tani (1974) allowed for the 
wake influence in connection with our instability problem) . 

For the base flow there is a consensus that the conceptual division 
into the "outer" irrotational flow and the viscous and rotational boun- 
dary layer and wake is justified on physical grounds even though careful 
mathematical treatment of the high-Reynolds number wake and its coupling 
with the upstream potential flow remains elusive. It is worth pointing 
out parenthetically a significant contrast between the irrotational base 
flow and our perturbation flow where vorticity (linearizable or nonlinear) 
is present upstream of the stagnation region. This complicates the matching 
between the approach flow and the perturbed flow near the stagnation line. 

ft 

This inner perturbed flow was conceived by Gortler (1955) and explored in 

j| 

much detail by Hammer lin (1955) , It is functionally equivalent to the 
Hiemenz component of the base flow and will be called the GH flow. 

The word "matching" above is used here in the sense of Van Dyke 
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(1964, p. 89) , i.e. applying systematic limiting processes which may 
assure a smooth continuation of the piecemeal approximations to the 
solutions of the Navier-Stokes eguaticxis. Such dovetailing between the 
Hiemenz flew and the surrounding potential flow has been generally assumed 
as possible but has apparently been tried more rigorously only by J. Cole 
(1968, p. 159) . Cole used the inviscid velocity potential and stream 
functicai for coordinates, but we fall back cn rectangular coordinates x, 
y, z, with z along the stagnation line, in order to illustrate several 
aspects of the matching process. But first we should recall some pro- 
perties of Hiemenz’ solution. 

b. The Hiemenz Scaling and Embedding 

Figure 1 (after Schlichting 1968) depicts some velocity distributions, 
streamlines and the boundary- layer thickness of the Hiemenz solution (plus 
the stretching of a material line B^^C^ into B^C^ for later purposes) . 

The figure also conveys the crucial fact that in this flew fields which ex- 
cerids to infinity in the x and y directions^ there is no characteristic 
length L and no characteristic velocity \! in terms of which distances x 
and y and velocities u and v could be measured or made dimensionless. This 
scaling feature leads to self-similar families of solutiens in terms of 
Hiemenz magnification variable X\=/a/V y, where a is an arbitrary velocity 
gradient (or an arbitrary ratio \^/i) and V the kinematic viscosity. The 
solutions (tabulated and graphed as functions of T] in most texts on viscous 
flows, e.g. on pp 89-90 Schlichting 1968)^ can be scaled to fit the inner- 
most field near the stagnation line of any irrotational flew approaching 
a blunt cylindrical body by proper choices of a and V. As 0 grows to 
infinity the viscous Heimenz solutions approach the potential flow 

Uroo = ax, Vr^ = -a(y - 6 ) , 6 = 0.648/ v/a , (1) 

•k 

where 6 is a constant displacement thickness due to the viscous decelera- 
tion near the wall at y = o, see Fig. 1. The asymptotic approach to the 
potential flew is very rapid so that already at ri = 2,4 the H flow is within 
1% of (1) , the definition of 5 in Fig. 1. The vorticity layer is thin indeed: 
the vorticity generated at the wall spreads by weak molecular diffusion and 
is convected back by the stronger mean approach flow. We should anticipate 
similar physics for our perturbed flow field. 

Figure 2 indicates the manner in which three inviscid potential flows 
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Cexpressed in complex notation with Z = x + iy') are related to the 
limiting Hiemenz flow CD • It is clear that in small neighborhoods of 
the origin, as Z-*-o, all three flows merge into CD if the x axis is 
shifted by the displacement thickness Cwhich vanishes in the limit v+o) : 
y* = y - 0 , For the cylinder, a = 2V^/R, and for the flat plate strip, 
a = V^/L, vdiile the wakes are disregarded. Their presence would primarily 
decrease the magnitude of the scaling potential velocity gradient at the 
the origin, a = unsolved problem being to 

determine by how much. CThe wakes also introduce u -velocity fluctuations 
at the stagnation line, dominated by the overall instability of the wakes. 
This larger-scale motion is not believed to influence the instability of 
the approach flow, although it is definitely present in most measurements 
such as in Fig. 9 , discussed in Section 6e.) 

Incidentally, the scaling parameter a also appears in the solution 
for the technologically important heat transfer rates in the stagnation 
region. Therefore, to provide a reliable reference base for discussion and 
evaluation of experimentally determined heat transfer rates, OU/3x)^ 
must be evaluated carefully. 

c. Outer Scales and Matching of Base and Perturbed Flows 

From the preceding illustrations it should be clear that the limit 
n -► <» does not correspond to free stream conditions. Yet this is an 
implicit assumption in the various mathematical manipulations of the 
instability problem as will be seen later. In our illustrations the 
variables Y y/R or y/L will be called the outer variable in terms of 
which the boundary conditions at the upstream infinity can be imposed 
and satisfied; Y -> «». The connection to the viscous boundary conditions 
at the cylinder comes through the adjustable self-similar H solution and 
its matching with the outer solution, see Cole (1968) p. 159. In simplest 
Van Dyke terms the matching requires that the limit of the inner solution 
as n coincides with the outer solution as Y 0. Here the inner 
solution would consist of the Hiemenz solution near the stagnation point 
and of the local boundary layer solutions for other azimuthal positions 
(where the outer variable would be a curvilinear coordinate equivalent of 

y 

/R) . Matching involves functional overlap over a region and the explicit 
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prescription Y -»• 0 emphasizes that this region is right at the wall . 

For the base flow with irrotational behavior upstream of the body the 
limiting process y -► presents no difficulties. Desirable applica- 
tion of the technique to reattaching shear flows with mean vorticity 
Cinger 1976, Ginoux 1958, 1971), however, could not well disregard viscous 
diffusion in the outer flow, cind the equivalent of the process Y would 
require extra finesse. An even more serious situation faces us in our 
instability problem where the up-stream perturbation is vortical and 
is subject to the vorticity-modifying straining field induced by the 
body. In the physical problems of interest free-stream turbulence decays 
steadily in the uniform flow before entering the straining, partly ampli- 
fying approach flow. The turbulence characterization as a fixed percentage 

of therefore has little meaning as a boundary condition for Y ->■ «». And 
if we seek true eigenfunction growth, somehow elicited by the presence 
of free-stream turbulence, as discussed in Appendix 1, where should we 
prescribe the vanishing of the perturbations? The condition of vanishing 
disturbances in the Hiemenz limit appears inadequate unless all 

elements of the instability mechanism oculd be shown convincingly to 
remain confined to the region right at the wall. The eigenvariable h for 
this viscous-conditioned problem is along the stream direction in contrast 
to the usual instability problems where r| lies across streamlines, see 
Appendix le. Wilson and Gladwell (1978) identify the instability as 
centrifugal and not only accept the p matching, but require that the 
perturbations vary exponentially in p, see Appendix le for details. These 
and other instability issues are discussed further in Section 7 after the 
relevant steady and unsteady experimental information is reviewed. At 
this stage, however, it is clear that the behavior of disturbances in the 
region corresponding to r) ^ « and to finite Y will require special atten- 
tion . 

d. The Hiemenz Flow Field 

Perturbation equations inevitably involve various characteristics 

of the base flow, in particular of the inner Hiemenz flow, as coefficients. 
For reference purposes we therefore recapitulate succinctly the salient 

features of the Hiemenz solution. 

Under the inner n = y magnification all two-dimensional stag- 
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nation flows with an irrotational outer stream ttirn out to be similar, 
irrespective of the shape of the body, as in Figure 1. For irrotational 
outer flows the stagnation streamline comes in perpendicularly to the 
wall, and under the magnification, it is difficult to distinguish between 
the wall and its tangent plane. The body shape (in particular the 
leading-edge radius) determines the strain rates at the origin, labeled 
a in Equation (1) , Introduction of a dimensionless stream function 
with appropriate scaling of the dependent variables: 

= ax (n) , Vjj = Cri) ;r? = y (2) 

automatically satisfies the continuity equation. The governing equation, 
that of the x mcmentum component, is expressible in terms of r) alone, for 
all outer irrotational flow fields^ e.g, Schlichting (1968) p. 89: 

+ $$••)+ 1 = 0 . (3) 

a n H n 

Here the highest order derivative describes the viscous effect as usual, 

the nonlinear terms (in parenthesis) reflect the inertial acceleration 

and the last term corresponds to the constant pressure gradient p/ ^ x 
2 

=j>a , made dimensionless. The no-slip boundary conditions at the wall are 

simply $' (0) = ^>„(0) = 0. The blending of the viscosity-conditioned 

n H 

region with the irrotional flow of Equation (1) is accomplished simply 
by the remaining boundary condition: 

1 

- — ( 0 ) •*“ 1/ as n (4) 

Because the Hiemenz scaling, Equation ( 2 ) , absorbed the usual free non- 
dimensional parameter Re^ characterizing the ratio of viscous and inertial 
terms, self-similar solutions became possible. As a result the corres- 
ponding instability equations and boundary conditions contain no base- flow 
characteristic parameter either, a parameter which could become "critical" 
as instability sets in, just as Reynolds number becomes critical in Toll- 
mien-Schlichting instability. 

For a circular cylinder of radius R, one may hopefully try to switch 
to cylindrical coordinates and consistently retain only the lowest order 
of an expansion in powers of 0 so that x = R0 and r = R(l+ j/v/a n/R) • If/ 
following Tani (1974) , one keeps the notation of U and V for the corres- 
ponding azimuthal and radial velocity components, one needs only to modify 
the second of definitions (2) to 
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and obtain the governing equation - which is identical to Equation (3) , 

No new R-based parameter appears to regularize the stability considera- 
tions, the solution $ (n) remains the same, and only the radial velocity 

” -1 
component is slightly modified by the factor (1+ y/R) for 0<y<y , 

m 

corresponding to 

The solutions for $ , $*, and are tabulated and graphed in 

il n H 

Schlichting C1968) , pp. 90-91. One should wish to supplement these with 
the solution of the y— momentum equation which yields the pressure field: 

P„(x,ri) - P„(0,0)-l/2p{ax^+vf (n) }-pVa$* . (6) 

U II £1 li 

The main pressure rise in the stagnation region is inertial, as expected, 

the last term in Equation (6) , due to viscous resistance, being very small. 

* 

The Hiemenz displacement thickness is 6 = 0.648 ^v/a and momentum thickness 

is 0.292fJ7^, 
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3. FREE-STREAM DISTURBANCES AND STEADY UPSTREAM VORTEX FORMATION 

3a. Spatial and Temporal Dependance of Free-Stream Disturbances 
The notion of instability implies small perturbations from established 
states of flow, i.e, from the matched approach- and Hiemenz flows of 
Section 2 in the present problem. There are an infinite number of possible 
perturbation fields; they can be usefully idealized into several characteristic 
types. Since there is no evidence that small acoustic disturbances do 
affect significantly the flow around cylindrical bodies (except at separation 
of the shear layers which initiate the wake) irrotational disturbances are 
not usually considered. At a given section y= const, of the approach flow 
the rotational departures from uniform flow are ftinctions of x and z 
which can also be (a) time-independent, as for Wcikes at very low Reynolds 

I • 

ntmbers, (b) time-periodic, as in the case of non-turbulent von Karman 
wakes, (c) random-stationary, as free-stream turbulence is usually portrayed 
in textbooks or (d) a superposition of motions (a) , (b) and (c) . 

Appendix If discusses the spotty and sporadic nature of decaying free- 
stream turbulence , which at least for the purpose of some instability 
modeling may perhaps be idealized in terms of separate, energy-sharing 
coherent structures. 

For small disturbances we also have Fourier decompositions at our 
disposal to reduce the complexity of the functional representation in 
X, z, and time. Thus for a low- Re wake , diffusing from an x - oriented 
upstream cylinder (i.e. perpendicular to the z orientation of our stagnation 
line and parallel to the line in Fig, 1 ) the free-stream perturbation 

would be AV = /(z) , resembling an error function, see the right side 
of Fig. 3b. The error function is expressible in a rapidly convergent 
Fourier series so that an analysis with 

AV = A C04 (2ttz/X) = A eO-4 k z (7) 

would be meaningful. Alternately, a single row of thin parallel wire s with 
mesh M = X, called a zither for short, can be shown to generate a decaying 
wake in the form (7) with A(y) decreasing exponentially over a significant 
distance, see Kellogg (1965) , 

These are of course but simplest illustrations of general Fourier- 
integral decompositions of arbitrary fields in the wave number space, which 
are singled out because of their special usefulness for our problem. The 
disturbance vorticity associated with (7) is parallel to the x axis and 
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has the magnitude -Ak 6'in fcz, A number of visualization techniques, 
e.g. Fig. 12 , have disclosed sporadic vortex formations in the stagnation 
region with a doxainemt content, as in Fig, 3a, even when the upstream 
disturbances were not generated steadily and preferentially as for Fig, 3, 
but by general turbulence grids at high Reynolds numbers. Thus Equation 
(7) seems to be on the right track and our main task is to quantify the 
mechanisms which select this orientation of disturbemces and which 
(more critically) prefer a relatively narrow reinge of wave numbers k when 
incoming turbulence contains broad spectra of frequencies and wave 
numbers . 


3b, Knowledge and Specification of Free-Stream Conditions 

II 

A quarter-century after Gortler devised the analytical formulation 
of the problem, cuid focused on disturbemces of the type (7) , the second 
part of the task has not been convincingly carried out and the 1978 voices 
of Wilson and Gladwell question the accomplishments in the first part. 

One reason for this is that until a few years ago the experimentalists did 
not seek the simplest controlled disturbemces of the type (7) , Instead 
they tackled directly the most complex unsteady problem of free-stream 
turbulence and tried to interpret it in terms of the steady formulation 
(7) . Without exception, the wave-number characteristics of the turbulence 
remained unknown and the complicated function /(x,z,t;k^,k^,k^,) of the 
preceding Section was specified by a single number, Tu = A\/) in 

percent. 

Miyazaki and Sparrow (1977) comment that poor specification of Tu 
(e.g. its y location and decay rate) and inaccuracy of Tu measurements are 
probably the most serious flaws which face the interpreters of heat-transfer 
enhancement. Furthermore, when Tu on the order of 3% or higher is claimed, 
nonhomogeneity in x and z of both Tu and can be suspected; the problem 
then also acquires nonlinear features from the beginning. Even for linearizable 
Tu values, the unsteadiness and the broad wave-number spectra would make it 
essentially impossible to ascertain whether part of the response represents 
a true free eigenfunction behavior, even if space-time correlating tech- 
niques were used, see Appendix lb and Ic. 

3c. Steady Vortex Formations Upstream of Obstacles 

Any student of vorticity and three-dimensional flows is well advised 
to ponder the lessons of the beautiful Frontispiece in **Inccmpressible 
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Aerodynamics", a remarkable collective effort edited by B, Thwaites (1960). 

The flow in the Frontispiece featuring a steady five-vortex formation 
upstream of a cylinder within a laminar boundary layer, is described 
as "radically different" from its two-dimensional counterparts and its 
"complexities" are intended as a challenge for the reader, A correspond- 
ing seven-vortex formation was studied by Norman (1972) by a single 
smoke-sheet generator which he moved across the boundary layer to trace 
the individual streamlines in the plame of symmetry as they wind their 
ways into the different vortices, e.g. in Fig, 4, There, two near-by streamlines 
with high total pressure are also outlined by dashed lines. Furthermore 
arrows are added to emphasize the jet-like character (evident in Norman's 
movie) of the upstream directed motion along two of the streamlines as 
they push against the decelerated fluid near the wall. 

In 1973, guided by this experience, the present author conjectured 
that steady separated vortices upstream of cylindrical obstacles could also 
foirni away from sidewalls, as in Fig. 3, if the streamlines with the slowest 
velocities of an oncoming steady laminar wake had sufficiently low stagnation 
pressure P^. Such upstream separation must entail a free stagnation point 

(Fig. 3) where the jet moving upstream from the stagnation point S on the 

cylinder is brought to rest at a total pressure P^ by the oncoming low- 

inertia fluid in the center of the wake. The stagnation pressure at the 

cylinder point S opposite the center of the wake must be lower than that 

farther along the stagnation line S^, where the full undisturbed P^ is 

achieved. There is therefore acceleration along toward S In addition 

to the usual Hiemenz acceleration "around the obstacle", a , in the x 

H 

direction perpendicular to the plane of the lower Figure 3. Apart from the 

viscosity effects, the magnitude of acceleration a tends to limit the 

H 

maximum velocity that can be generated along the stagnation line toward 

S, and therefore P^, The pressure at the free stagnation point, P^, 

should be very nearly p^+ - ^V) ^ , if p^ and represent "initial" 

conditions two or three diameters upstream of the cylinder (often denoted 

by subscript i) . For conditions under which p builds up beyond P-, a 

s t 

reverse flow and vortex formations should ensue as in Figure 3 , 

3d. Steady - Flow Experiments of Hodson and Nagib 
By studing wakes of single euid multiple cylinders at low 
Reynolds niunbers as these wakes approached either rectangular or circular 
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cylinders, Hodson and Nagib C1975) verified the present author's conjecture 
of Section 3c and provided approximate criteria for the rather sudden appearance 
of the vortices such as the one seen in Fig. 5a. They also clarified considerably 
the effect of free-stream unsteadiness on the vortex behavior and on heat 
transfer rates at various locations of a large cylinder, but these observations 
will be discussed later. In Fig. 5a, the dye marks the center of a steady 
subcritical horizontal Wcike which impinges upon the vertical face of am 
elongated rectangular cylinder. By madcing the length of the main cylinder six 
times its width, 0=23.75 mm, the upstream influence of its wake (end of Section 
2b) was minimized so that the bifurcation of the oncoming wake at (Fig. 3) 
and its tight roll-up into a pair of horseshoe vortices was steady to the eye. 

In the side view of Fig. 5a the horizontal traces of the two vortex cores obscure 
the fact that the dye did not penetrate the boundary layer of the target 
cylinder although some of the fluid from the edges of the wake should do so. 

The flow was always observed simultaneously in the plan view through a 45° 
mirror arrangement. The reader is referred to Hodson and Nagib (1975) for 
discussion of the plan-view observations which revealed significant alterations 
of separation ^rid wake behavior of the main cylinder in the proximity of the 
vortices. While the implications of such observations for various flow 
characteristics, including heat transfer, are of interest, for the present 
purposes the important findings are the existence of steady vortical response to 
steady upstream nonhomogeneity and what Hodson and Nagib call the threshold of 
incipient formation of the vortices . 

3e. Evidence of Steady-State Threshold and Implications for Theory 

Figure 6 presents an ad hoc empirical correlation of the conditions which 
separate states with and without vortical cells in several series Of visuali- 
zation experiments by Hodson and Roadman (1975) under Nagib' s guidance. 

They were done with smoke in air and dye in water (in a water channel and a 
towing tank) and with single and multiple wakes which impinged on one 
rectangular and one circular cylinder. The experiments were carried out with 
the concept of the vortices as a module in heat-transfer convection rather than 
with stability theory in mind. 

Other attempted correlations involving wake momentum defect and wake 
vorticity met with lesser success. One should keep in mind that these 
quantities and the thickness 6 of the stagnation boundary layer 
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eu:e very difficult to measure without probe interference at the conditions 
needed for visualization and therefore had to be estimated. In Figure 
6 , 6 /d = cCv/DV^)* 5 , c = 6.59 cind 4.66, gives the estimate for the boundary 
layer thickness at the stagnation line of the rectangular and circular 
cylinders, respectively; M is the "mesh" or lateral distance between wake- 
generating multiple cylinders, which provide a variation of Av without 
change in the important lateral length scale after the individual wakes 
merge; Cf is the grid solidity d/M; and the subscript of Re indicates 
the length on which the Reynolds number was based. The ordinate (Ay/d) 
is related to the strength of the centerline momentum defect empirically - 

Hodson and Nagib (1975) , Nishioka and Sato (1974) , Kellogg (1965) . The 

★ 

term C<5/d) in the abscissa is proportional to D/6. One would wish 
that the correlation be more physical but even in this steady laminar 
case too many parameters enter the formation mechanism. Without theoretical 
backing, there is no reason then to expect the threshold state to map 
into a single curve in a two-dimensional representation. 

How does this threshold as perceived by Hodson, Nagib, and Roadman 
relate to instability theory in the light of the discussion of nonlinear 
effects and diagnostics in Appendix Ih and li? Let us assume that the 
system would allow initially exponential growth in time as its characteristic 
Reynolds niomber passes the neutral condition, if the disturbances 

could be kept infinitesimal. This assiimption is made in most of the 
theories thus far. Sutera et al (1963) , Sutera (1965) , Williams (1968) 
and Inger (1977) also start with infinitesimal disturbances but suppress 
time dependence in their analysis. Let us also assume that free-response 
finite disturbances in the system could be characterized by a single 
amplitude A(t) and that the system has am inverted bifurcation at so 

that its neutral curve has a shape as in Sketch 1. These assumptions 
allow for the simplest threshold behavior in temporal stability. It is 
the behavior anticipated by Tani (1974) on the basis of the heat-transfer 
variations with intensity of free-stream turbulence of Kayalar (1969, Fig. 1) . 
When the free response is characterized by A and a single scale . say a 
dimensionless wavelength X*, a X* aixis out of the paper has to be added 
and neutral states then correspond to surfaces in three dimensions rather 
than simple curves in Sketch 1. The more basic parameters there are, the 
higher the dimensionality of the neutral states. 

* Since 6 is proportional to >^Dv/U^. 
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According to bifurcation theory of the free eigenfunction behavior^ 

e.g. Joseph (1974) , the dashed part of the neutral curve separates the 

domain of stable initial system states (Re,A(0))at the left from the 

domain to the right, where the states evolve unstably. The unstably 

developing states to the right of B-Re^^ may reach the stable amplitudes 

of the upper part of the neutral curve B-C provided turbulence does not 

set in - there is disagreement on the stability of states BC to all 

disturbances, in temporal instability, any initial state (Re, A(0) ) to 

the left of B“Re decays until it settles on the stable line 0-Re of 
cr cr 

zero disturbance energy. The system follows a path such as 1', if Re 
remains constant. 

Actually there were two distinct experimental procedures: Hodson's 

data in Fig. 6 with a horizontal slash through the symbols denote experi- 
ments in which the wake-casting cylinder was fixed, at y=yj^/ while was 
slowly increased, and his data with a vertical bar through the symbols 
signify that remained constant while y^ increased until the vortices 
disappeared. According to the simplified theory, as increases in the 
first set, the system would either decay along paths like 1', or, if A(0) 
were on the B-Re^^ dashed line or to the right of it, the path would 
shoot up to the stable equilibrium B-C, as along path 1 in Sketch 1. Thus, 
if the temporal nonlinear growth rates past the neutral curve B-Re^^ were 
fast enough, all the experiments of the first set in the region marked 
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"CELLS" in Fig, 6 should be observed to map into the curve B-C in Sketch 1, 
with amplitude A dependent only on Re and not at all on the initial ampli- 
tude . Corresponding to the region marked "NO CELLS" in Fig, 6, the 
experimenters should see not just the absence of cells but rapid temporal 
decay of all normal mode disturbances to the zero energy states corres- 
ponding to Sketch 1, To change experimentally entails 

inevitable time lags. Therefore, identification of the initial time t 

o 

for the system (Appendix lb) is not feasible for the experiments and one 
should not expect to make definite conclusions from whatever transient 
temporal behavior might be observed. 

The experiments unquestionably indicate the presence of nonvanishing 
disturbances thus contra-indicating behavior 1’ in Sketch 1; in fact the 
steady disturbances even grow spatially over part of the region , most 
likely in accordcuice with mechanisms discussed in Sections 5 a and 5b. 

Could it be argued that these nonvanishing, spatially growing velocity 
variations represent a forced rather than a free response of the system 
to spatial initial conditions of Hodson et al, within a temporal framework ? 
Forced responses are not dealt with in the insteibility theory and are not 
present in Sketch 1. But if the interpretation in terms of forced motions 
were correct how could one maintain that the nonlinear threshold behavior 
which results from strengthening of this forcing is transformed into a 
free response describable within the framework of Sketch 1? Since the 
computable free eigenfunctions of the linear theory would be modified in 
an unknown manner by nonlinearity in their evolution toward the time- 
independent modes corresponding to the neutral curve segment B-C, there 
seems to be little hope of operationally confirming the free-eigenfunction 
interpretation in this case. Esthetically , the forced-response interpre- 
tation would appear self-consistent. Given future computational capabilities 
it might even be numerically verifiable starting with a time - independent 
nonlinear perturbation of the Hiemenz flow with constant finite A values 
in the upstream boundary condition, Equation (7) . The aforementioned 
papers of Sutera et al (1963) , Sutera (1965) and Williams (1968) in fact 
proposed such a forced, steady, weakly nonlinear model and carried out 
some rather limited computations of uncertain convergence. At present, 
neither these computations nor the experimental information are adequate 
to validate the forced response viewpoint in any detail. 
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Thus far the comparison of experimeaits in set 1 with behavior 
portrayed in Sketch 1 dealt with paths of the type 1 ' , We have already 
seen that corresponding to paths of type 1 the observable amplitude A 
should depend only on the Reynolds number. No assessment of the strength 
of the nonlinear vortices , such as the positions of the free stagnation 
point in Fig. 3b or 5a, were made in the experiments of type 1. 

However experiments of type 2 , where Re was kept constant and the forcing 
amplitude gradually decreased by moving the wake-generating rods upstream, 
do partly bear on the same issue of exclusive Re-dependence of the equilibrium 
amplitude A past the threshold. To conform to the nonlinear behavior of 
free eigenfunctions in Sketch 1 , the amplitude A^ would have to remain 
unchanged on B-C as the rod position y^ increased until a state at or 
past the unstable neutral curve segment B-Re^^ was reached - at which 
time the amplitude should snap from A^ to a finite value and thoi decay 
along the path 2 in a manner simular to 1'. (The snap from A^ to below 
B-Re^^ would allow for the possibility of hysteresis in this peculiar 
open system.) 

Again, no specific record was made of amplitude characteristics as 

y^ was increased^ but in type-2 experiments amplitudes (as characterized 

by the position of S^) are directly and immediately comparable. According 

to the memories of the experimenters, A^ did not remain constant during 

procedure 2 , again contra-indicating behavior 2 and further undermining 

the basis of the two-dimensional Sketch 1. However, as y, of a single 

1 

rod varies, not only the wake disturbance amplitude, but also its scale 
varies (while for contiguous multiple wakes X* remains constant): . In a 
three-dimensional generalization of Sketch 1 to A, Re, X* space, 
procedure of type-2 for a single wake would correspond to a neutral curve 
in a plane Re=constant starting with A^ and varying with X*. The 
experimental observations do not therefore rule out a theoretical threshold 
behavior in a generalized three-dimensional space. 

According to the author's simplified reasoning in Section 3c, based 
on upstream forcing plus enhancement due to adverse pressure gradient 
(Section 4a) , one would expect slow growth of the vortices (after they 
are fully formed) , as long as the effective strength of the disturbances 
can be increased. These expectations are not inconsistent with the 
Hodson-Nagib-Roadman experiments. The general impressions seem to be that 
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some forcing response is present, and that somehow the characteristics of the 
boundary layer on the cylinder enter the quantitative determination of the 
geometry of the horseshoe vortices . Whether these features cem be 
incorporated in some modifications of the present linear euid nonlinear 
theories will be teiken up in Section 7. 

Figure 3b shows a fully formed double vortex which, according to 
Hodson, corresponds to conditions "sufficiently away from the correlation 
curve" in Fig. 6. Hodson also comments that "usually a small hairlip 
is observed at the threshold" (illustrated in Fig. 8 of Nagib - Hodson 
1977) . This corresponds to a thin cusp of dye extending upstream from 
S in Fig. 3b when two dye lines are used on either side of the oncoming 
wake*. If this visualization discloses faithfully all the upstream 
motion, there would exist a second type of free stagnation point (within 
the boundary layer?) besides the regular type in Fig. 3b, albeit over 
a very limited range of parameters. Also, developments past the threshold 
would not be discontinuous but rather rapid as a function of the steady 
disturbance strength. 

In accordance with Gen. Eeshotko 4 one would like to see more detailed 
independent corroboration of the findings of the Hodson-Nagib experiments, 
especially since they were not devised to test stability theory. As stated 
in the Introduction, new information also gives rise to new unresolved 
issues; the formulation of these issues in this Section should help to focus 
the new experiments. 


* It is possible that the slight momentum defect from these two far-upstreeim 
slender ejectors was sufficient to modify the critical conditions at the 
^®^t)~momentum stagnation line and caused the cusp. When the dye source 
is in the wake-generating rod as in Fig. 5a, the above interference is 
absent. 
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4 DIGRESSION ON OTHER INSTABILITIES LEADING TO STREAMWISE VORTICES 


(4a) SiTnilarities and Differences between Three Flow Families 
It is instructive to compare the growth of the vortices discussed in 
Section 3e to that of Taylor vortices in an annulus between two cylinders 
when the inner cylinder rotates, e.g. Stuart (1971, p 348) , and to the 
growth of Gortler vortices in a boundary layer with concave curvature, e.g. 
Bippes (1972) , Tani (1962) , and Wortmann (1964) , Usually, the stagnation 
flow is linked to the Taylor and Gortler centrifugal instability , e.g. 

Stuart (1963, p 505) , because the streamlines in Fig. 1 also have concave 
curvature. We shall see that the analogy may be weak at best (Section 5d) 
and that there may be significant differences in the stability character- 
istics of the three flows. Already the geometries of the base flows differ 
substantially: there is no "upstream", "downstream", or "free stream" in 

the annulus, while all these stream limits are present in the boundary layer 
on concave walls, and "free stream" and "upstream" coincide in the stagnation 
flow. One might expect the structure of the eigenvalue problems to reflect 
the associated physical differences, but the original formulations in terms 
of temporally evolving normal modes were essentially identical in the 
three cases. In particular, the growing vortices were all expressed as 
afj( 9 ) exp(3rsin(k_5) , alternating in the spanwise direction J5 • Also, the 
eigenvariable was always taken directed across the viscous layers of the 
base flows. However, in the region of validity of the stagnation Hiemenz 
flow this 9 direction (Fig 1) is essentially upstream rather than cross- 
stream, a key physical difference which is hidden by the formal similarity 
of the eigenvalue problem, 

(4b) Couette Flow Instabilities 

The temporal formation seems suited to the physics of the Taylor vor- 
tices in the closed annulus domain . Any local vorticity disturbance is not 
convected away by the stream (as in the concave boundary layer or in the 
stagnation flow) but remains within the closed confined space. Linear 

theory for an annulus of infinite length indicates that past a critical 

2 

Taylor number, T (T=(V.d/v) »(d/R.); R.,V., radius and velocity of inner 
cr 1 111 

cylinder; d, gap between cylinders) normal perturbation modes of computable 
wavelength grow as exp/3? . A number of weakly nonlinear models (Stuart 

1971) suggest that the higher terms are stabilizing so that a steady vortical 
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motion, with A*^A / should result. The centrifugally stimulated motion 
cr 

feeds repeatedly back on itself and the strong constraint of mass conserva- 
tion leads to the topology of an infinite series of vortex cells . 

Experimentally, when the inner cylinder of an annulus of finite 
length L is suddenly started rotating at speeds somewhat above vor- 

tices spread axially inward from the ends, grow very rapidly past their 
linearized stage, and indeed reach nonlinear ly limited finite amplitudes. 
However, the pattern evolving after the spin-up assumes a non-unique span- 
wise wave length fairly near that of the theoretical linear for L-^oo. 

Its exact value evidently depends on the initial distribution of vorticity 
in the domain at the time of the sudden spin-up (see Appendix lb for com- 
ments on initial value problems in closed regions). Snyder (1969, p 288) 
concluded on the basis of imaginative experiments with initial disturbances 
that "the mode having the highest amplitude initially (always) prevails". 

While the rapid initial growth in time appears consistent with the 

linearized exponential growth (if one keeps in mind that the impulsive 

acceleration is finite and that instriiments have resolution thresholds) , 

the vortex spacings A/d after two dozen revolutions will exhibit some 20% 

nonuniformity , presumably related to the initial distribution. Snyder 

(1969, Fig 6) showed that asymptotic equalization of the individual vortices 

depends on slow spanwise diffusion through the full chain of cells, with 

2 

acceptable adjustment time of approximately L /6v'. (Unless use is made of 
high-viscosity fluids "steady state is not reached for periods of from hours 
to days", and therefore "all past measurements of A are open to criticism", 
Snyder, loc, cit. , p 286) . Snyder's findings were essentially confirmed by 
the overlapping investigations of Burkhalter and Koschmieder (1974) and 
Cole (1974) . 

When the other experimental procedure of slowly creeping up to is 

employed, the initial velocity perturbations consist of (a) the flows induced 

at the axial boundaries which do not move in accordance with the ideal 

Couette velocity distribution, V=Ar + B/r, expected far away from the ends (but 

see Coles and Van Atta 1966) , and (b) the nonuniformities induced by the 

final speed increment of the inner cylinder. The latter influence is 

probably small, but the boundaries bring forth clearly discernible "anomalous" 

vortex pairs at the ends at about 0,8V and continue to induce additional 

cr 
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"ghost" vortices as is approached. However, for sufficiently long 
annuli, "the speed at which all divisions between cells just become visible 
has been found to agree well with the calculated critical speed," Cole 
(1976, p 7) , Cole (1974, 1976) and Burkhalter-Koschmieder (1973) provide 
us with detailed systematic observations and error estimates. The L- inde- 
pendence of the final rapid growth of the center cells in long annuli 
probably justifies viewing the final vortex formation process as a true 
instability which may grow nearly exponentially when the center disturb- 
ances are still small, in accordance with Taylor's linear normal-mode 
theory . 

The vortical foirmations induced by the ends are part of the initial 
conditions for the last ^V, increment, but since the no-slip conditions at 
these boxindaries persist, the sidewall motion may be considered as consti- 
tuting steady forcing by the boundaries . This sidewall action thus bears 
functional similarity to the role of the upstream rods in the Hodson-Nagib 
experiments. (There is even scattered evidence that at very high speeds 
the first seeding of randomness and turbulence starts in the forced, anoma- 
lously larger, vorticity formations at the ends, in the Couette case, just 
as it unquestionably starts in the forced wakes in the Hodson-Nagib and 
other experiments.) Thus, in presence of the steady forcing and the 
initial conditions (e.g., the initial flow distributions of Snyder), modes 
of motion similar if not identical to Taylor's linear normal modes are 
elicited away from the ends, grow and settle to nonlinearly modified stable 
modes with wavelengths near the unforced theoretical linear A 

cr 

As grows beyond an increasingly broader family of normal modes 

is admissible according to the linear theory. However, a small subset, 
judged by measured A values, is actually observed. Nonlinear theory for 
infinite L (Stuart, 1971, p 354; Stuart and DiPrima, 1978) suggests that 
some linearly admissible modes are nonlinearly unstable and that resonance 
mechanisms among sideband modes and the harmonic of the fundamental are 
responsible for the narrower range of allowable A values. For finite L the 
additional requirement that the nxmaber of cells has to be an integer, imposes 
strong quantization restrictions on the observed modes and leads to hyster- 
esis effects in experiments with increasing and decreasing V^. These major 
features of nonlinear Taylor flow behavior probably hold a lesson for anyone 
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expecting a well-defined A variation for the stagnation-flow vortices. 

In line with oxir interpretation of the end-wall geometry as forcing, 
Koschmieder (1975) introduced 0-rings, 0.17d in height, at equal distances 
along his annulus in order to impose and stimulate preferential growth 
of vortices on the flow. And stimulated they were: "the first photograph 

with signs of vortices along the entire colijinn was taken at 0.03T^^", i.e. 
at 0.17V^^ ccanpared to approximately 0.8V^^ for the flow with normal end 
effects (Cole 1976) . In quasisteady procedures, vortices with 
1.47A^^^ ^ ^ unstable and split up into two vortices in a non- 

critical, continuous manner. Consistently with nonlinear theory, the range 
of stable vortices 0.69 ^1.18 (roughly the same for quasisteady and 
sudden-start procedures) was narrower than the admissible range of linear 
modes; see Pig 2 of Koschmieder (1975) . In accordance with Gen. Reshotko 4 
(section 1), one would wish for corroboration of the nonlinear theory with 
nonpermanent forcing elements, e.g. vortex generators,, which could be with- 
drawn during the experiment once V had been exceeded. It should be the 

cr 

nature of this closed nonlinear system, as it is the implication of Snyder's 
experiments (1969), that the forced components of those periodic motions with 
within the allowed set would decay after the withdrawal of the generators, 
i.e. the motion would evolve into the allowed A^ nonlinear mode without 
rearrangement. Those forced motions with A^ outside the allowed set would 
break up and rearrange themselves ultimately into one of the allowable set. 

In contrast , according to the Hodson-Nagib experiments, it is the nature of 
our stagnation- flow field that removal of the forcing rods upstream of the 
cylinder removes the horseshoe vortices near the stagnation line altogether . 
The disturbances simply are convected out of the open system and so are 
their response formations . 

In terms of bifurcation theory the configuration for the amplitude A (t) 
of the equilibrium disturbance states of the system with infinite L is super- 
critical , as in Sketch 2; see Stuart (1971) or Stuart (1977) . Below T^^, 
the quiescent states denoted by the heavy line equilibriiam 

states and other perturbed states should decay in time along paths like 

path 1' in Sketch 2. Past T the quiescent states along the dotted hori- 

cr 

zontal axis correspond to unstable equilibrium. Staijle equilibrium states 
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Sketch 2 


are now associated with the regularly bifurcating branch with defi- 

nite non-zero amplitudes. Other pure initial states would approach equili- 
brium along paths like 1, Again, Sketch 2 describes states completely 
characterizable by a single amplitude A(t;T);and to accommodate specifica- 
tion of wave lengths \ or wavenumbers k we must generalize the equilibrium 

curve T -B to an equilibrium surface in the A,T,A or the A,T,k space, 
cr 

The dots past B in Sketch 2 allow for higher bifurcations , specifi- 
cally for the so-called wavy instabilities {Stuart 1971, Fig 2) at higher 
Taylor numbers (which depend sensitively on the ratio of the cylinder 
radii) . The annular vortices develop an integral number of waves in the 
azimuth which propagate around the annulus so that a fixed probe records a 
discrete frequency (and weak harmonics) . As an example. Coles (1965) 

obsejTved 3 to 7 waves for R./R of 0.874 and a length- to-gap ratio of 27.9. 

1 o 

As these waves set in^ they presumably first grow exponentially in time until 
they reach new nonlinear equilibria where the vorticity distribution is 
characterized by terms like exp [i (me-i2TrftJ in addition to the earlier 
factor sin(k3+ const) . We shall see that Hassler's 1972 experiments on 
the instabilities in the stagnation region also indicate a time-periodic 
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behavior of the streamwise vortices which has not been previously discussed 
as a possibility » Section 6e, 

The bifurcation diagram of Sketch 2 strictly corresponds to the infinite 
annulus and does not allow for forced responses, steady or unsteady. Never- 
theless, according to our brief account, forced motions are a fact of life 
in any physical realization of the flow. The concept seems to provide 
fruitful guidance for experimental diagnostics even though the forced motion 
cannot be generally treated as a small perturbation. In bifurcation theory, 
attempts to accommodate finite geometry (but not other forced motions) are 
sometimes described as studies of imperfect bifiircations . According to p 141 
of Stuart's 1977 review, "it is believed (though in no sense proven) that a 
smooth transition takes place" and that the new equilibrium states correspond 
to a curve in Sketch 2 which runs somewhat above, but follows the heavy equili- 
brium lines rather closely, with the maximum deviation in the neighborhood of 
T^^ . Stuart goes on for three pages to describe the difficulties of match- 
ing this abstract theory to experiments and the existing discordant views of 
both , 

(4c) Gortler Vortices in Concave Boundary Layers 

As pointed out in Appendix Id, the temporal formulation of the growth 
of streamwise Gortler vortices in concave boundary layers definitely does not 
fit the experimental observations: there is no time variation in the early 

development and the growth is in x . To the author's knowledge the determina- 
tion whether the instability is subcritical or supercritical has not been 
made , theoretically or experimentally. We therefore do not know whether to 
anticipate behavior associated with Sketch 1 or Sketch 2. In fact, there 
have even been svibstantial discrepancies in the computed neut ral curves for 

If 

the critical Gortler number G versus wave mmiber k (G = Re if ; R, wall 

t n ' R 

curvature; th, momentum or displacement ^ickness or l/vx/u , depending on the 
author) . A survey of the problem and a collection of contending neutral 
curves is found in Herbert (1976) . Herbert pinpointed inconsistencies in the 
difficult treatment of streamline curvature as a major contributor to the 
discrepancies. He also concluded that the gradual thickening of the boundary 
layer should be a first-order effect, of more importance than for the TS 
instability. An extended abstract of Floryan and Saric (1979),, available at 
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this writing, suggests that they succeeded in devising a systematic non- 
singular approach to the curvature problem via the method of multiple scales 
and the DiPrima-Stuart (1972) scaling of disturbances in streamwise vor- 
tices. They also traced some of the past discrepancies to computational 
errors. Interpretations of experimental results should thus probably start 
from the new curves of Floryan and Saric, 

Experimental study of the Gortler vortices is itself very difficult 
because measxirements of steady perturbations in space are much more insensi- 
tive and inaccurate than those of unsteady perturbations at a fixed position . 
Nevertheless, the hot-wire studies led by Tani (Tani 1962, Tani and Aihara 
1969) leave no doubt that there is at first an exponential growth in x which 
diminishes downstream either due to nonlinear saturation or due to detuning 
of the normal modes in a thickening boundary layer. Wortmann's (1964) 
ingenious Telluriimi visualization technique in water generally corroborates 
the findings of Tani and coworkers. He also observes that the normal and 
spanwise velocity perturbations are twenty to thirty times smaller than that 
in the streamwise direction, thus providing support for the DiPrima-Stuart 
(1972) scaling. The rate of rotation of the streamwise Gortler vortices is 
apparently exceedingly slow (completing barely a 90® turn before transition) 
in comparison with our stagnation -flow vortices, such as in Figs 5 and 12 , 
also in water. Stereoscopic studies of hydrogen-bubble traces in water by 
Bippes (1972) indicate that the slow rotation ultimately accumulates bulky 
ridges of low- inertia fluid along directions where the vortex flow is away 
from the wall and that a complicated, nearly time-periodic secondary 
instability develops in proximity of these ridges . For comparison of its 
characteristics with features of secondary instabilities of nonlinear TS 
waves and of the wavy-form instabilities of Taylor vortices, see Bippes' 

1972 discussion, pp 165-175. 

There is some disagreement among the experimenters concerning the 

nonlinear developments and the nature of the breakdown (which are iinlikely to 

be universal, anyway) . However, there is general consensus about the linear 

behavior. At incompressible speeds without heat transfer, the amplification 

rates are relatively slow* in comparison to TS amplifications. As the 

* Gortler instability is suspected of contributing to transition on the con- 
cave walls of supersonic wind tunnels . 
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experiment of Bippes in Appendix la indicates, the more highly amplified 
normal modes are evidently selected from disturbances present in the layer. 
However, the maximimi of these cimplif ication rates with respect to wave nijinber 
may be shallow. Thus by driving a disturbance with a wave length twice the 
preferred one, by winglet devices above the boundary layer, Tani and Aihara 
(1969) completely preempted the normally observed mode. Presumably, the 
forced mode was built up to nonlinear levels before the usual mode was of 
any consequence and thereafter became nonlinearly inhibited, "Mild selecti- 
vity to wave length" would explain why the wavelengths observed in foxir 
different experiments of Tani and Aihara were "almost independent of free- 
stream velocity". Evidently, small nonuniformities of flow inherent in a 
given experimental arrangement can force the Gortler response in their image 
provided the wavelength is amplified by the instability mechanism. 

(4d) Reprise of Comparison of Instabilities with Streamwise Vortices 

In Section 4a we listed a nxamber of a priori similarities and dissimi- 
larities of the three families of flows with streamwise vortex formations, 
which in textbooks are linked by susceptibility to centrifugal instability. 
All three were originally foimulated as temporal instabilities. The experi- 
mental evidence reviewed in the preceding sections strongly suggests that 
only the Taylor instability in Couette flow can be considered temporal in 

II 

nature and that the Gortler instability is definitely spatial. Even though 

end-effects cause end vortices to be formed below the linearized critical 

conditions in Taylor instability, at the vortices all along the annulus 

grow rapidly — to undoiibtedly nonlinear configurations. Spanwise equalization 

however, takes place rather slowly by diffusive processes in an adjustment 

2 

time on the order of L /6v , In the other instabilities such times are not 
available as the vorticity is transported downstream. For these, we should 
therefore expect (and indeed do see; Tani 1962) spanwise variations in ampli 
tude and wavelength of streamwise vorticity unless the vorticity is uniformly 
seeded upstream. 

Among the three flows, only stagnation region vortices exhibit velocity 
components opposite to the direction of the momentum which drives the flows . 
It is in the development of this local counterflow that the threshold-like 
phenomenon appears. Neither Taylor nor Gortler vortex behavior seem to pro- 
vide any guidance for understanding of the Hodson-Nagib experiments and in 
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particular of the threshold-like behavior in Fig 6. Nor is the model of 
subcritical instabilities in Sketch 1 truly consistent with the observations, 
incomplete as they are. 

However, in all three systems various aspects of the observed linear 
and nonlinear behavior seem to require for adequate explanation the presence 
not only of instability modes but also of forced response to initial and/or 
complete boundary conditions which are present in any given physical realiza- 
tion, This is the message of the Taylor-flow experiments of Snyder (1969) , 
Koschmieder (1975) , Cole (1976) , Biirkhalter and Koschmieder (1973) and of the 
selective stimulation of Gortler vortices with upstream grid and with narrow 
streamwise heating elements in the wall by Bippes (1972) , and with winglets 
placed outside of the boundary layer by Wortmann (1964) and Tani and Aihara 
(1969) , Quantitative mathematical elucidation of the details of the mecha- 
nisms constitutes the generalized problem of linear and nonlinear receptivity. 
Some of the forced responses may decay in time or space, after determining 
the early development of the normal mode response, as in Bippes' case of the 
small heating elements, while in other cases the forced response remains a 
factor throughout. Mathematical clarification of the characteristics of 
various categories of possible behavior would be helpful . The nonlinear 
aspects seem to be rather prohibitive at present, but are needed at least 
conceptually. With hindsight of the preceding sections, we can conjecture 
that the steady Hodson-Nagib vortices represent in part a forced nonlinear 
response which may or may not involve a free response, i,e. non-zero solutions 
of the associated homogeneous system of equations and boundary conditions. 

The disappearance of the vortices as the wake-producing rods are moved 
farther upstream seems to support the conjecture. 
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5 3JECHANISMS AND LINEARI2ED DISTURBANCE EQUATIONS 
5a Prandtl*s Inviscid Linearized Formula 

That there should be a tendency to accentuate steady^flow defects 
between neighboring streamlines^ of Fig 3^ in regions of adverse pres- 

sure gradients has, of course, been Icnown for a long time. That it could lead 
to local reversal of the flow in absence of a ’’sidewall^’^ with its no-slip con- 
dition (which was present in the case of Fig 4), that is intriguing. Prandtl 
(1932) considered small steady flow velocity differentials and V 2 between 
two stations where stream velocities are and and the stream tube 

areas A^ and A^ , Assuming that the process of acceleration or decelera- 
tion was rapid enough to neglect the slow action of viscosity, he used the 
conservation of stagnation pressures and of mass flux to derive a first order 
expression for the relative magnitudes of the defects; see Corrsin (1963) ,p. 553 , 



Here it is the change in the stream tube areas which reflects the role of the 
given streamwise static pressure drop or rise, p^- p^ . Actually, Equation 
(8) has been used primarily for assessment of channel contraction ratios A^/A^ 
which would depress steady-flow nonuniformities to a tolerable level . It is 
equally valid for expansion of stream tubes in front of cylinders which takes 
place in planes perpendicular to their axes , The stream divergence accord- 
ing to the asymptotic Hiemenz solution for the stagnation region is depicted 
in Fig 1. With Az in the direction of the cylinder axis, a typical stream 
tiibe area expands to A^ = Az and Equation (8) can be 

rewritten as 


LV 


2 



^ 2^2 


AV 


1 * 


(9) 


While the result (9) basically rests on the fact that a given increase 
in pressure along streamlines has much more effect on fluid elements with 
lower inertia, it can also be interpreted in terms of the companion variable, 
vorticity, because AV^/Az -- and A^ 2 /Az -- 60^^ • Equation (9) 

then represents the increase in vorticity by stretching of vortex lines, 


which is the mechanism usually associated with the increased velocity 


29 



differences in stagnation regions. However, Equation (8) is clearly valid for 
non-planar geometries, say upstream of asymmetric three-dimensional bodies, 
where the corresponding analysis of the vector velocity field would be unwieldy, 

5b Forcing by Diffusive Vorticity Convected from Upstream 

Figure 7, borrowed from Kestin and Wood (1970) , presents schematically 
computational results of a steady second-order theory as well as the coordinate 
system, matching that of Fig 1, i.e, with z axis along the stagnation line 
and the free stream coming from y ^ + oo . Weak periodic steady laminar wakes 
-Vj^cos kz (such as generated by Hodson and Nagib's two-dimensional grid in 
Fig 3) and the free stream, -V^ , convect vorticity = -kV^sin kz from 
upstream "infinity" toward the cylinder which is approximated locally by its 
tangent plane y = 0 in Fig 7 . 

As long as the stream remains unaffected by the pressure field of the 
cylinder, roughly for y ^ 2.5D , the vorticity perturbation decays exponen- 
tially by diffusion - exp{“[vk^ (y^- y)/VQ] } from some initial reference 

location y^ , For uniform base flow, constant, the sources in the vor- 

ticity equation (A2-5) of Appendix 2 are all zero and the vorticity evolution 
is governed purely by diffusion. This viscous diffusion and decay of vorticity 
does not cease, of coiirse, as the cylinder is approached. It complicates the 
formulation of the upstream boundary conditions for the disturbances. The 
formula suggests that the disturbance amplitudes may have to be prescribed at 
specific values y^ on the order of 2D if the free-stream forcing is to be 
taken into account realistically. This is of course way beyond the reach of 
Hiemenz variable ^ . 

There is a degree of confusion concerning the outer behavior in this 
direct inflow problem. To the experimentalist the subsidence of vorticity 
with decreasing y represents a most natural decay , To stability analysts 
looking outward from the body, the behavior appears as an exponentially 
divergent growth of disturbances to infinity. Such a divergence is not con- 
sidered permissible for the normal modes which are presumably rooted in the 
body boundary layer. In the Hiemenz plane, where the ve locity of the base 
flow is allowed to diverge proportionally to rj as r} "O , the diffusive 
behavior of the vorticity is further camouflaged by straining. Asymptoti- 
cally algebraic decay (or growth) then replaces the exponential decay obtained 
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for uniform - 

As discussed in Appendix the normal modes do not develop out of 
nothing, but in conjunction with the forced response of the system, The outer 
diffusive behavior of the vorticity represents predominantly the forced res- 
ponse, For the Hodson*>Nagib experiments, analyzed at length in Section 3d 
and 3e, the relation between the forced and free responses remained elusive, 
but the necessity of paying attention to the forcing was clear , Figure 7 , 
kindly supplied by Professor Kestin, shows explicitly the upstream boundary 
conditions (which were deleted in the published version) and thereby under- 
scores the peculiar nature of the instability in the direct inflow stagnation 
region. What looks like a non-homogeneous forced problem is reducible to a 
homogeneous normal -mode problem by switching the boundary condition 

= constant to the condition on the derivative, dV^/dy = 0 as \Tf ^ oo * 
Assuming for the moment that the full problem was solvable in the r) variable 
alone, the full solution would presimiably consist of the forced solution plus 
a dovetailing solution of the homogeneous equations with = 0 as ^ r 
if such a solution exists. The mechanisms underlying the equations (which 
differ only in one boundary condition) would in this case be the same . 

5c Vorticity Diffusion and Vorticity Stretching 

In approaching the cylinder the flow diverges in the yx plane perpen- 
dicular to the cylinder axis and the convected 6iJ^ vortex lines are stretched 
by the Hiemenz flow as already seen in Fig 1. For instance, the vortex line 
segment ^ 1*^1 stretches into segment convected with the 

main flow and thereby generates a vorticity source in the sense of the right 
hand side of Equation (A2-5) , According to Equation (A2-5) , the rate of 
change of observed as one moves with the material particles toward the 

cylinder depends therefore to the first order on the race between the dominant 
stretching source oj (9U/3x) (equal to -oj (3V/9y) by continuity) and the 
diffusion loss rate v(3^W^/3y^+ d^U)^/dz^) , The diffusion due to the cvirva- 

ture of the vorticity profile in the yz plane is governed by the square 

2 2 2 2 2 2 

of its wavelength X : V(9 oj /9x ) = -Vk Oi^ = -v(4tt /X )0) , The squared 

length controlling the diffusion v(9 ti)^/3y'^) is that associated with the 

approach to the cylinder, i,e, usually a much larger term, on the order of 
2 

D , except in the immediate vicinity of the cylinder, roughly up to twice the 
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height of the vortices in Fig 7. Furthermore, above this height this loss, 

2 

besides being small, is of the same sign as "Vk A slightly conserva- 

tive estimate of the ratio of the competing local stretching and diffusive 
rates is therefore -X /v4tt or +X /v4tt every\rfiere in a narrow 

wedge-like domain centered on the yz plane npstream of the body, where the 
above approximations are valid. 

Far upstream, ^V/ ^y is negligible, the above ratio approaches 
zero, and we recover the exponential decay. The gradient ^V/^y scales 
with ^0^^ since total is itself and the only characteristic 

length is D, For a circular cylinder in potential flow, ^V/^y grows 
to approximately 1,2 Vq/D at y = 0,25D and reaches its maximum of 
4V^/D at the edge of the boundary layer. There, however, the diffusion 
term V ( caJ y^) cannot be neglected so that the maximum VIT^D 

X 0 

of the approximate ratio of the stretching rate to the diffusion rate may be 

misleading. In particular, we do not expect the ratio to increase indefinitely 

with wavelength A . As A increases, the x-dif fusion decreases, so that it 

will cease to be a controlling factor. Furthermore, as A increases, the 

likelihood of building up the stagnation pressiire P beyond P , the 

5 f 

free stagnation-point pressure, also decreases — see the discussion in connec- 
tion with Fig 3b and the dual view at the end of Section 3c. 

Clearly, in case of countertrends, general physical considerations such 
as in the preceding paragraphs need to be replaced by actual solutions of the 
governing equations. As will be seen in Section 7, such solutions have not 
yet yielded convincing criteria for the "most amplified wavelength A ", 

The preceding discussion at least brought out the importance of both the 


diffusion and vortex stretching as the dominant physical factors in the steady 
phenomena portrayed in Fig 5a. It also suggests that we can expect net 
vorticity growth, solely from vortex stretching, especially very near the body. 


5d Centrifugal Destabilization 

Despite references to concave streamlines, no comparable demonstration 
of the role of centrifugal destabilization has been reported in the literature. 
According to Equation (A2-5) , there are no additional terms in the Gortler- 

II 

Hammerlin vorticity equation beyond those identified above, not even within 
the boundary layer. Wilson and Gladwell (1978) conclude on the basis of a 
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generalization of Rayleigh’s inviscid instability criteria that "the 
(centrifugal) destabilizing forces are confined to the region where (the 

base flow vorticity) is appreciable, namely the boundary layer". The same 
conclusion can be drawn from Stuart's remark (1963, p, 505): "In such a flow, 

there are certainly regions where the circulation (product of local velocity 
and local radius of curvatxire) increases as the local center of curvature is 
approached normal to the curved streamlines; consequently Rayleigh's 
criterion suggests instability," In accordance with Helmholtz' theorem this 
change of the circulation vanishes outside the boundary layer. In other words, 
curved irrotational flows are neutrally stable. 

Wilson and Gladwell identify two types of terms in the x and y 
momentum equations as due to centrifugal acceleration. However, they conclude 

U H 

that these are relatively smaller and end up with the Gortler-Haininerlin system 
of equations . 


u 

5e The Gortler Perturbation Mode, Its Temporal Dependence 

In 1955, long before any visual evidence of the formation of horseshoe 

II 

vortex pairs at blunt two-dimensional leading edges was available, Gortler 
proposed that unstable, vortical, spanwise-periodic perturbations might account 
for the reported increase in unsteadiness and heat transfer in that region. 

The linearizably small three-dimensional perturbations were superposed on the 
Hiemenz flow, Equations (2) , so that the velocity and pressure fields were 
represented essentially as follows: 


U = ax{$y(ri) + (n)e^^ cos k?} 

W = -^ W^(n)e^^ sin kC 

P = Pj^(x,Ti)+ pvaP^(ri)e^ cos k^ , 


\ 

> ( 11 ) 




Here and henceforth T = at and k stands for the dimensionless wave- 
number related to the dimensional spanwise wave length X by k=(2'iT/X) /v/a 
since the spanwise coordinate z undergoes the same magnification as did y, 
namely z = >/^/a ^ . Substitution of (11) in Navier-Stokes equations and 
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linearization yields a system of jaiomentum equations and the continuity 
equation 


Undifferentiated P 2 (u) appears in the z momentum equation and P^(T)) in 

the y momentum equation. Differentiating the former and eliminating Pj^(r^ ) 
leads to a sixth order differential system in the velocity components which 
is expected to be subject to six boundary conditions. Mathematical proper- 
ties of the system and its solutions form a subject beyond the present dis- 
cussion, much as they influence the question about their relation to the 
physics of the problem. However, the simplest, nearly self-contained ecfuation 
is the vorticity equation, which is well worth examining here. 

We first note that Gortler allowed for non-propagating temporal insta - 
bility through the sole unsteady factor exp ( /3 7- ) which cancels out in the 
continuity equation. The eigen-number ^ enters all the other equations; in 

particular in the vorticity equation below it enters as f3 CaJ and acts as a 

2 

sink of vorticity, always in unison wxth the z - diffusxon sxnk k if 

the equation is interpreted in a steady system. As discussed in Section 3, 
the experimental evidence of Hodson and Nagib does not support the exponential 
time dependence, but neither does it definitely exclude its possibility. 
Furthermore, hot-wire evidence of Hassler , Figs 9 and 10 and Section 6e, dis- 
closes quasl-periodlc time variations in response to random turbulence Tu 
of approximately 0,85% measured at y^^ = 0.83D upstream of the cylinder 
(10 cm or 120 mesh-length downstream of the screen) . Such experimental evi- 
dence contains unknown contributions from forced motions , At any rate , the 

II 

Gortler normal-mode system does not allow for such temporal behavior. 

5f The Vorticity Equation and the Domain of its Validity 

Its mathematical properties and physical interpretation have been cen- 

II 

tral to most analyses, starting with Hammerlin (1955) , who studied in detail a 
function g( *7 ) without identifying it as vorticity. If we enlarge the 
definitions (11) by = aco^(ri)exp(3T) sin , its relation to V-W velocity 
components becomes ^^1 through the continuity equation (12) 

to as well (V^- k^V^- U^)/k , Hammerlin ’s g function. Since the 
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perturbations are assumed to be small one can obtain an equation for 
directly from Equation (A2-5) in Appendix 2: 

“i ■■ <*'h" “ ° • 

The vortlcity stretching source rate u)^(3U/3x) and the diffusion loss 

rates and V(3^a>^/3y^) , used for the less restricted approach 

to the cylinder in Section 5c, become for this specific perturbation pattern 

2 

(which is immersed in Hiemenz flow) * ^®spectively . 

(The term represents the observer derivative as he follows the mean 

Ulemenz path toward the cylinder.) The ratio of sources to sinks becomes 
» independent of when k^ is not too small and rt > 6 for y 

diffusion to be negligible. In contrast to section 5c this lower bound on a 
characterization of the net local amplification increases monotonically as the 
cylinder is approached. 

According to Section 5c, the free-streara disturbances first decay due 
to vorticity diffusion before vorticity stretching reverses the trend; (see 
also the experiments of Sadeh, Sutera and Maeder , 1970b) , The specification 
of the boundary conditions for the system and its consequences will be taken 
up in Section 5g; however we already saw in Section 5b that the problem of 
computing the response to the prescribed spanwise-periodic inhomogeneity at 
"infinity" of Fig 7 is included in the homogeneous problem when the deriva- 
tive of is made to vanish at infinity. The monotonic rise in the net 

local amplification associated with the vorticity equation (13) indicates that 
the system cannot duplicate the occurrence of the minimiam where vorticity 
diffusion is balanced by vorticity amplification. This again reflects the 

ti 

fact that the domain of validity of the Hiemenz and the Gortler-Hammerlin 
equations should be measured in small multiples of the boundary layer thickness 
rather than in terms of the body characteristic length D and that application 
of the method of matched asymptotic expansions, MMAE, should be considered; 
see also Sections 2a - 2c. 

It is not surprising, therefore, that the study of asymptotic solutions 
of Equation (13) for tj ^ ao (which is not the same as for y -»■ , as 

brought out in MMAE) , remained the center of theoretical investigations as well 
as of contention, e.g, Kestin and Wood (1970 ) f Tani (1974 ) f lida (1978a, b) . 
Sadeh, Sutera and Maeder (1970a) in fact abandoned attempts to characterize the 
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approach to free stream through the magnified variable and used "cylindrical 

2 2 

coordinates" x,y Ci,e. 8 or :x /y very small as in our discussion lead- 
ing to Equation C5) ) , to obtain a Reynolds number parameter in the equiva- 
lent of Equation (13) for the outer y-field. However, no rigorous matching 
with the inner field as characterized by Equation (13) was attempted, Sadeh 
et al (1970a) oversimplified their system so that their Equations (22) and 
(23) are inconsistent and their subsequent results at best qualitative. In 
view of the sensitivity of these issues, it is worthwhile to review the 
vmtranslated Gortler-Hammerlin results and the discussion of the outer limits 
in their variable , 

5g The Outer Boundary Condition on 

The no-slip boundary conditions at the wall U^(0) = V^(0) = W^(0) = 0 
present no difficulties but the outer limit *7 oo does , on two covuits . As 
previously noted, the first is philosophical; do we view the overall problem 
as an open system where small but nonzero upstream pert\arbations on V of 
the form exp ( ) . cos (with ^ possibly zero) can be admitted in or 

beyond the limit >7 ->■ cO ? This would be the view of experimenters and of 
engineers concerned with observed enhancement of heat transfer. Or do we 
ask that all perturbations strictly vanish as 17 oc and then arise 
"spontaneously" along the paths of the fluid elements? This would be the 
pure normal approach confined to the Hiemenz scales. In Tollmien-Schlichting 

If 

waves and Gortler vortices, the normal-mode amplitudes are linked to the 
amplitude of the disturbance environment, though the mathematics of the 
linkage remains obscure. If a similar linkage existed here, solutions with 
both types of outer behavior would presumably be present. 

The second difficulty with the outer limit is associated with the 
singular nature of the 7 limit in which the base flow velocity itself grows 
without bound: as ^ 00 , Since this limit is fictitious (unless 

MMAE is applied) , what is a logical physical requirement for perturbations 
around it? 

Gortler (1955) first asked that all perturbations U^, V^, and 

approach zero as -v co , but then observed that the milder requirement 

V, /V^ — ► 0 should be sufficient, even if V, should increase indefinitely, 

10 1 
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5h Hammerlin'a Results and Some Implications 

The solvability of Gortler's system of equations and boundary condi- 
tions was investigated by Hammerlin (1955) , It involved first an assess- 
ment of the number of non^divergent linearly independent asymptotic solutions 

of the asymptotic equations in U and V in which S was replaced by 

IX H 

Y = - 0,648 and by unity as Y and <jO , Hammerlin then 

examined the feasibility of matching smoothly at some the outer solu- 

tions (made up of the available, linearly independent outer asymptotic 
solutions) to the well-behaved inner solutions, 19 - (made up of six 

regular linearly independent solutions) and the possibility of simultaneously 
satisfying the wall boundary conditions. For the continuous wave number 

2 ff 

spectrum 0 < /3 + k cT 1 , Hammerlin found an extra degree of freedom for 
the coefficients of the linearly independent solutions. Equivalently, for a 
chosen positive value of /3 + k < 1 , one "amplitude" can be chosen arbitrar- 
ily at , and then there would be a unique set of solutions such as those 

of Inger (1974) in Fig 17, For (^ 4 k^) > 1 , the extra degree of freedom 
disappeared because a previously convergent solution diverged and had to be 
discarded by setting its coefficient to zero. Hammerlin considered the 
possibility that the resulting nine homogeneous equations in the nine constants 
could have discrete eigensolutions but concluded that this could only be 
settled numerically. In a limited numerical search with H.J. Maehly on a 
1955-type computer no eigenvalues were found. 

For the so-called neutral, i.e. steady, solution, ^ = 0 , Hammerlin 's 

asymptotic solution for V involved algebraic variation with the dominant 

^ 2 

term ( 9 ” 0.648) . This would be allowable for k < 2 according to 

Gcjrtler's weaker boundary condition ^ as >9 ^ e>o , However, 

iJ 

according to Hammerlin the boundary condition V' 0 (needed to make 

^ 2 

^ 0 through Equation (12)) requires that k < 1 so that the above 

term would also decay to zero , Thus for perturbations of Hiemenz ' flow as 

K) CO , 0 implies -*• 0, According to Inger (1974 and 

private communication) , that is the case for his solutions such as those in 
Fig 17, 

The implication of this is that the forcing boundary condition of 
Fig 7 cannot be satisfied with ^ 0 within the Hiemenz flow, where 
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itself diverges as C *9 0,648) , Wilson and Gladwell C1978) believe that all 

instability jnechanisms should occur within the Hiemenz scale and reject 
Hainmerlin^s algebraic approximation to as unphysical. They maintain that 
any linearizable vorticity influence upstream from the boundary layer must 
decay exponentially outward with , As discussed in Section 5b, this 
algebraic solution has the earmark of diffusive vorticity strained in the 
convective Hiemenz motion downstream toward the body. If so, discussion in 
Section 5f and Hammerlin’s conclusion concerning 0 as >7 oo 

indicate that responses to unavoidable oncoming disturbances in this open 
system and their possible linkage with free normal modes cannot be accom- 
plished in the Hiemenz plane and that at least a two-region asymptotic match- 
ing process may be required to improve the correspondence of Gortler's model 
with experimental observations. Even then nonlinearity and unsteadiness 
(quasi-periodic or randomlike) are likely to stand between the theoretical 
model and the physical phenomena. In this writer's opinion, however, the 
main stumbling block is the confinement to the Hiemenz plane. It is for this 
reason that he does not discuss the painstaking^ weakly nonlinear treatment of 
lida (1978a, 1978b) , It does not seem to improve the physical picture, 
presiamably because it is locked in the variable. 
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6 EXPERIMENTS WITH UNSTEADY DISTURBANCES 
6a Possible Mechanisms of Enhanced Heat Transfer at the Stagnation 
Line 

As noted in Section 3b, all the reported experiments other than that 
of Hodson and Nagib had unsteady free-stream perturbations which clearly do 
not correspond to the assumed perturbations of Equation (11) . Most investi- 
gations, largely motivated by experimental observations of heat transfer 
rates in excess of the theoretical (laminar) predictions, Schuh (1953) , 
naturally focused on the technologically important cases at higher Reynolds 
nvimbers where the oncoming stream is usually txirbulent. However, the exis- 
tence of steady vortical flow fields, such as in Figs 3 and 5a, suggests 
that increased heat transfer rates could be due to (a) increased "mixing" in 
steady nonlinear horseshoe vortex formations associated with steady nonhomo- 
geneity of the oncoming flow; (b) increased heat transfer rates across tempor- 
ally thinner shear layers associated with some unsteady aspects of the flow; 
or (c) additive or more than additive effects of simultaneous conditions (a) 
and (b) . Small perturbation theory of Fourier-decomposed unsteady stagna- 
tion boundary layers yields only small second-order heat transfer effects, 
e.g. Lighthill (1954), so that effect (b) alone is unlikely to account for 
the observations. 

6b Visual Evidence of Hod son and Nagi b and of Sadeh and Coworkers 
It is instructive to watch the movies which Hodson and Nagib made while 
gradually introducing unsteadiness in the oncoming low-Reynolds niamber wakes 
such as in Fig 5b. One finds that even though the modulation of the vortex 
formations is relatively slow, the relaxation times in this through-flow 
system are extremely short in comparison with those in Taylor vortices; see 
Section 3g. The sequence of frames shown in Fig 8 was recorded at the inter- 
vals of time shown in the figure. Here the steady wake of Fig 5a was replaced 
by a periodic but still laminar von Karman vortex street at Re^ of 90 by 
increasing the speed threefold. At 111 diameters downstream of the wake 
generating rod the fluctuations are rather weak, e.g. Kovasznay (1949) , and 
this small, relatively high-frequency vorticity variation appears to change 
the mean stagnation region vortex cells very little. (For lower forced 
frequencies the effect could well be more substantial.) 
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Special interest centers on the response of the "amplifying system” to 
a moderate local disturbance marked P whose progress through the system is 
documented in Fig 8 and is described by Hodson and Nagib as follows : "The 

vortex flow module is always present on the average but is kicked out of 
position into the 'free stream' by the perturbation. When this occurs a 
large amount of the fluid in the vortices, which is marked by the dye, is 
ejected away, providing a very efficient mechanism for the mixing of the flow 
in the stagnation region with that in the free stream. It is conjectured 
that this mechanism also leads to added enhancement of heat transfer from the 
bluff body." The break up of the vortex pair and the ejection are seen in 
stages 5 and 6 of Fig 8. They are rapid and nonlinear and should give rise 
to some very thin temporal boundary layers. In Fig 8, the "disturbance" P 
marks a relatively small change in the oncoming unsteady but rather regular 
vorticity pattern so that the events can be clearly identified, contorted as 
the streaklines are. Apparently, when the disturbances are strong enough to 
cause mildly unsteady horseshoe vortex pairs to form we can expect sharper 
vortex interactions as the oncoming vortex pattern changes . Hodson and 
Nagib also mention interaction with vortices in neighboring cells when 
multiple wakes are observed. 

Sadeh, Brauer, and Garrison (1977) and Sadeh and Brauer (1978) present 
visualization studies using a single smoke filament injected just 0.32D 
upstream of the stagnation line. The short length of the history-integrating 
smoke streakline facilitates somewhat the interpretation of the motions in 
the region of interest which is viewed from the side, from the top, and from 
the rear. However, the momentum defect from the injection tiibe and its 
support system (enhanced by the adverse pressure gradient) is convected 
directly at the stagnation line and interferes with the sensitive motions 
over half of the span. The Reynolds number Re^ of the wake generating rods 
was 650, yielding a much more random onflow than that of Fig 8. Since the 
rods were only 18,3 mesh-lengths upstream of the main cylinder, there was 
ample upstream nonhomogeneity (category a) and unsteadiness (category b) , 

Each study describes detailed analysis of 62 frames of motion picture film, 
each frame with elaborate interpretative sketches . A companion motion 
picture, available on loan from NASA ^ makes one appreciate the complexity of 
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the erratic response of the rolled—up vortex structure to oncoming turbulent 
eddies and the possible interaction with the neighboring horseshoe vortices. 

At these upstream wake conditions visual information concerning the oncoming 
eddies (such as that in Pig 8, where the disturbance P is identifiable) is 
not obtainable, which hampers the interpretation of the responses and inter- 
actions , 

This writer has some reservations concerning a few passages of less 
critical interpretation* and the smoke ejector interference. The reader would 
benefit from studying the reports first-hand and could then judge for himself. 
Even if the reservations were correct, the net evidence amply demonstrates 
that increased turbulence intensifies distortions and vortex interactions 
such as described above in connection with Fig 8. In particular, in the top 
views, (Fig 11 of Sadeh and Brauer corresponding to sequential views of Y 
streaklines in our Fig 5b) , "an additional feat\ire of utmost significance 
hitherto unexplained is unfolded". To this observer, the sequence of 
stills and the associated sketches in their Fig 11 seem to portray variants 
of strong interaction between the one visible half of the dancing horseshoe 
vortex pair, the invisible oncoming upstream-modulated vorticity and probably 
the invisible neighboring vortex pair. The strongest interaction does not 
appear dissimilar in nature to that in the last three frames of present 
Fig 8 (where, however, the neighboring vortices are altogether absent) . 

Whatever the best perception of the motions may be , their very comp 1 ex - 
ity r aises the question of adequacy of analytical modeling . These motions can 
hardly be considered quasisteady and mildly nonlinear. Could they possibly 
be portrayed with satisf actoiry accuracy by the first two or three teirnis of 
poorly convergent, weakly nonlinear approximations grafted on steady .neutral 
solutions (/3 = 0) of the GH equations (11)? Do we have to adopt statistical 


* Seldom do two independent observers agree fully on the interpretation of 
single^ history-integrating streaklines in randomly modulated flows with 
strong vorticity. Also some of the disagreements may be a matter of termi- 
nology. The one important issue, which is discussed in Appendix 3, concerns 
the criteria for calling the boundary layer on the large cylinder turbulent in 
the classical sense. 
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techniques from turbulence theory because of the randomness? Is the forma- 
tion of the horseshoe vortex structures at all e3q>lainable within such 
theories? 

6c Controlled Onset of Unsteadiness and Change in Heat Transfer 

While the steady nonlinear cellular motion pictured in Figs 5a and 7 
may bring forth heat transfer rates in excess of those in a regular stagna- 
tion boundary layer, the unsteady motions described in the preceding section 
are likely to enhance the heat transfer even more. In other words, the 
observed increased heat transfer rates in presence of free-stream turbulence 
appear to be more likely due to category (c) in Section 6a than to category 
(a) alone. Hodson and Nagib devised a heat transfer experiment which 
covered the change from steady to unsteady controlled wake disturbances. The 
control was accomplished by moving the wake generators from outside while the 
tunnel was running and thereby inducing or suppressing the nonlinear horse- 
shoe vortices at will. 

They constructed a two-dimensional zither-like mesh of fine wires 
(diameter of 0.127 mm or 0.005 in,; solidities 0.05 and 0.10) and cast their 
wakes on an instrimented circular hollow pyrex cylinder (diameter of 50.8 mm 
or 2 in.) in the arrangement of Fig 3b. As the wind-tunnel speed was slowly 
increased these wakes experienced transition to an unstable state revealed by 
an oscilloscope trace from a continuously monitoring hot-wire anemometer. 

This occurred near Re^ of 25, depending on the proximity of the zither wires 
to the adverse pressure gradient generated by the cylinder. A gold strip 
less than 10 mm thick was baked onto the front of the cylinder, covering 

C . 

angles -11 ^ ^ 11 . Whenever the hot wire would indicate instability of 

the wakes, the strip-chart recorder measuring the temperature of the gold 
strip would register a sharp decrease. Thus, substantially increased cooling 
by the air stream took place when unsteadiness set in. Since the vorticity 
released into the wakes of the "zither" wires increases continuously with 
Reynolds number, even in presence of an essentially two-dimensional instabi- 
lity, and since the strength of steady vortex cells, such as in Fig 5a ^also 
increases continuously with Reynolds number, the sudden change in cooling 
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can be ascribed* to the onset of unsteadiness . 

For methodology and detailed quantitative results the reader is 

referred to Hodson cind Nagib (1975) or Nagib and Hodson (1977) . Their 

results relevant to flow instabilities are as follows: for steady wakes of 

the zither, the average Nusselt number over the gold strip increased by 5% 

at Re^ of 19 and by 10% at Re^ of 24 when the zither was brought to the 

y/d range of formed vortex cells. As already mentioned, a discontinuous 

increase in Nusselt number occurred in this y/d range when multiple Karman 

vortex streets were sensed near Re, of 25. At Re, of 27. the increase of 

a d 

Nusselt number due to combined effects of the cells and the Karman-street 
unsteadiness was about 23%. At Re^ values of 31, 37 and 42 the rise leveled 
off at approximately 29%, with increased scatter. The experiment demonstrates , 
therefore, that the formation of steady vortex cells does increase heat trans - 
fer rates in the stagnation region, category (a) , and that presence of rather 
regular xinsteadiness increases them more , category (c) . We should keep in 

mind that the wakes at these low Re, values below 42 are not turbulent. The 

a 

conditions probably correspond closer to the unperturbed patterns of Fig 8 

than to the disturbed pattern after the perturbation P interacts with the 

vortex. These are difficult experiments and corroboration of the results and 

their extension to higher Re , values would be desirable , 

d 

6d Experiments in a Towing Tank 

Our account would not be complete without discussing and digesting the 

most detailed experimental information thus far on the vortex structure in a 

• ** 

turbulent environment, that from Gortler's Institute for Applied Mathematics 
and Mechanics in Freib\irg, Germany. The reader really should examine the 

4/ V 

German originals of Hassler (1971) and Colak-Antic (1971) to appreciate the 
foresight, special technique, and care it takes to elicit the small-'scale 
details of this unstable and vinsteady phenomenon in vivo. Its unsteadiness 

* Evidence from continuous changes in variables (velocity increase and de- 
crease past the critical condition; back and forth motion of the wire grid 
without stopping and opening the wind tunnel) is of course more convincing in 
diagnostics of instability effects than from discrete changes in variables 
(such as installations of different grids). Discrete changes in parameters, 
which require stopping of the wind tunnel, invariably bring about unmonitored 
changes in some conditions of the flow or instrumentation and increase 
scatter substantially. This and the need to fair discrete data impair omr 
ability to doctiment cause and effect in instability phenomena from discrete 
changes in variables . 
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dictated the use of hot-wire anemometry; the need for X magnitudes on the 
order 6-15 mm for adequate spatial resolution and visualization called for 

, V 

water as a medium. The hydrogen-bubble visualization experiment of Colak- 
Antic and the hot-wire explorations of Hassler were therefore carried out in a 
water tank 0,23 x 0,28 x 3 m. Water provided relatively high Reynolds numbers 
and visualizability at low speeds, 5-40 cm/s, but at the price of (a) 
special techniques and constant recalibrations because of the conductivity of 
the medivim and (b) the need to move the cylinder model with all instrvnnenta- 
tion (including turbulence screens, traversing mechanisms, the hydrogen- 
bubble cathode wire and the camera) on a sliding carriage along the channel. 

Two details are of particular relevance to the data to be examined : 

(1) the hot-wire signal itself was used as a guide to minimize the vibrations 
induced by the carriage motion so that these did not influence appreciably 
the 0-80 Hz spectrixn of interest; (2) the bare 90% platinum-rhodiimi hot- 
wire sensors had a diameter of 5 microns and lengths of 1 - 2 mm and were 
welded onto conical prongs of like material. In Figs 9-10 and 13-16, due to 
Hassler, the wavelength was 7 + cm (i.e. the wire averaged the signal over 
roughly one-fourth of a vortex cell), the boundary layer was 1,8 mm thick 
(Fig 14) , and the closest probe position was 0.5 mm from the cylinder, 12 cm 
in diameter. Thus when we consider the results in the plane of symmetry 
(from which the angle f is measured) , we should again allow for the possi- 
bility of some probe interference ; the flow is instability-prone and the 
probe does force upon it additional boundary conditions including an addi- 
tional momentimi defect. For larger ‘f x/R, the component U of the base 
flow (Fig 1) builds up and the danger of upstream interference rapidly 
decreases. 

Since the wire senses only the velocity normal to itself, the same 
Fig 1 or Equation (1) also helps us to assess the relative wire sensitivity 
to the U and V components of the velocity. The flatter the streamline, the 
larger is the response to U, a fact to remember in interpreting the figures 
because Hassler always calls the signal u, irrespective o f local stream line 
orientation . 

6e ^ Unexpected Quasiperiodic Behavior 

Figs 9 and 10 of Hassler disclose the presence of various types of 
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unsteadiness r which over the 0,16 second intervals (2,4 cm of cylinder 
travel) may display quasiregularity, some of it of unJcnown origin. These 
fluctuations are conditioned by free stream tiurbulence characterized by Tu 
of .85% (same for all the figures of Hassler) but unfortxmately not by a 
spectriim; see discussion in Section 3b, Even though the 12 cm cylinder had 
a 24 cm long splitter plate attached at its rear, the wake oscillated at 
approximately 8 Hz (see lowest of Figs 10, past separation) . For cylinders 
without splitter plate, upstream influence of both sets of vortices in the 
alternating wake produces appreciable fluctuation near the front stagnation 
line, at doiable the wake frequency. The small sample at 15.6 Hz in the 
upper trace of Fig 9 for ^ = 0 ® suggests that the effect may appear even in 
the presence of the splitter plate . 

Because of the acceleration around the cylinder, the boundary layer 
thickens very slowly so that little error is committed if one uses a single 
scale £ ~ 1.8 mm, even at ^ = 45°. Then the indicated y positions can be 
converted to the Hiemenz scale: rj'^1,3 - 1.4y, for rough comparison with 
Figs 17, 18, and 20 of eigenfunctions of the steady linear theory. The 
linearized vortex height corresponds to = 0. The trace at y = 5 mm in 
Fig 9 with the largest unsteadiness corresponds then to approximately twice 
that height where the steady perturbations have dropped 30 - 40% from 
their maximum. These traces exhibit a remarkably regular mysterious frequen- 
cy of about 86 Hz (60 - 70 Hz in other longer samples, according to Hassler 's 
private communication) . It is helpful to recall that Fig 9 portrays the 
behavior in the stagnation plane so that small vorticity perturbations are 
convected from the region of the lowest trace against the pressure gradient 
to the region of the uppermost trace. The lowest trace corresponds to t) 
of about 27, i.e. to the asymptotic region of the Hiemenz flow beyond the 
confines of Figs 1, 7, and 18. 

As we move away from the interference-sensitive stagnation plane to 
'f = 15°(x<v,1.6 cm~8.8^ ) the mysterious frequency is still present, randomly 
modulated, according to the first trace of Fig 10. Here the signal would 
represent mostly the temporal fluctuation in U, if it were comparable in 
magnitude to that in V. The remaining three traces make it clear that the 
temporal tangential fluctuations damp out farther around the cylinder: the 
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boundary layer remained laminar despite its spanwise corrugations cuid 
despite the relatively high free-stream turbulence of 0.85%. 

*rhe higher frequency oscillations of Figs 9 cind 10 are representative 
of the quasiperiodic behavior mentioned in Section 4b in connection with 
wavy instabilities in Couette Flow. The turbulence level of 0.85% here is 
undoubtedly much lower than the unspecified level of Sadeh and Brauer, who 
saw no quasiperiodic patterns. However the same level of turbulence is high 
enough to hide the periodic Tollmien-S chi ich ting response in non- stagnation 
boundary layers in most investigations, Hassler's diagnostics seem con- 
vincing enough to exclude as possible causes the relative vibrational 
motion between the hot-wire and the cylinder-attached velocity gradients, 
and purely electrical effects. Unsteadiness there should be as a response 
to the grid turbulence, but what mechanisms would selectively amplify a 
narrow band of frequencies just at the outer edge of the boundary layer? 

it 

The linear instability formulation of Gortler does not allow for such temporal 
periodicit y. Could it possibly be associated with secondary instability of 
the nonlinear vortex formations as in the wavy instability of Taylor vortices 

■r 

(Section 4b) and the Gortler vortices (Section 4c)? Unfortunately, the 
plans to verify that the effect is indeed generic and to study its parametric 
dependence in a larger towing tank at Freiburg have fallen through with the 
closing of the laboratory. 

V' 

6f Visible Patterns of Colak-Antic 

V V’ 

Colak-Antic 's two snapshots of the hydrogen biibble "time lines". Figs 
11 and 12, display the "top view", vertically downward onto the cylinder, 
with the leading edge seen tangentially at the upper edge of the photographs . 
The corresponding velocities of 11.9 and 19,4 cm/s bracket the 15 cm/s 
velocity of all the Hassler data in Figs 9-10 and 13-16, The lines of bxibbles 
are released at equal intervals at *f = 11** = 0.19 rad, y = 2 mm, i.e. just 
beyond the edge of the boundairy layer and away from the sensitive stagnation 
plane. Once formed upstream of this station, streamwise vortices should 
follow the pattern of Figs 7 and 20 even if they "dance around a little". 

If the bubbles were released at the height of the vortex centers where = O , 
a short time At later the bubble line Ay = V^^'At would resemble the profile 
of the curve in the lower Fig 20 except that the bubble density would 
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increase or decrease depending upon the local narrowing or widening of the 
stream tubes in th.e y — 2 crossplane. For bubble release above or below the 
height of the vortex centers the bubble time lines are also subject to span- 
wise distortions due to motion indicated in the lower Fig 20. Ultimately, 
had they had enough time, most bubbles would migrate toward the vortex cen- 
ters according to the streamlines of the upper Fig 20, At the speed of 
11.9 cm/s of Fig 11 the upwelling motion between the vortex pairs (near 
cos kz = 1 of Fig 20) is clearly visible and the ft shapes suggest that the 
bubble line was released below the vortex centers. The bubble lines also 
distort in the streamwise direction because of the variation of mean and 
perturbation tangential velocities U(y) as is especially visible in Fig 12. 

At the higher Reynolds number of Fig 12 the irregularities seem to be more 
pronounced, including a tendency for stronger vortex pairs to interact and 
even absorb neighboring weaker vortex pairs. Some of this behavior is not 
as surprising as it might bverwise be after the lessons of the Hodson-Nagib 
movie associated with Fig 8. 

6g On t he Causes of Spatially-Fixed Mean Patterns in Nominally 
Homogeneous Turbulence 

If free-stream turbulence were truly homogeneous one would expect 
vortex pairs to arise and disappear irregularly along the stagnation line as 
the different free-stream vorticity formations sweep in and out of the stag- 
nation region, giving instantaneous impressions similar to that of Fig 12. 

If the amplification along the paths had a sharp peak in wave number for the 
more complex unsteady perturbation equations here operative, one would expect 
a preferential wavelength to occur rather frequently but randomly along the 
span . However, true homogeneity would dictate no spanwise preferential 
location of vortex pairs, not even on a nonlinear basis. It is difficult to 
imagine an upstream influence from a finite-size vortex structiore (which is 
being washed out of the stagnation region) strong enough to force the inde- 
pendent oncoming free-stream vorticity formation to replicate its pattern. 

That oncoming vorticity formations are rather powerful even in presence of 
strong mean vorticity fields was illustrated in Hodson-Nagib, Fig 8, And yet 
substantial spanwise regularity in the mean does frequently occur in nominally 
hcanogeneous fields, as the remaining figures of Hassler will demonstrate. To 
the present writer this suggests that even subtle steady departiires from 
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spanwise homogeneity in the oncoming stream do get picked up in an amplify- 
ing process and enforce a steady mean pattern, (The effect could be likened 
to that of a minute mean drift making itself felt in face of random molecu- 
lar motion, except that in this analogy the drift would be amplified.) When 
occasionally a rare experimenter investigates the character of turbulence 
downstream of a turbulence grid or screen, with care great enough to resolve 
small mean changes in presence of large A.C. signals, he invariably finds 
that there is more inhomogeneity and nonisotropy than expected, e.g. Bradshaw 
(1965) and Grant and Nisbet (1957) . In particular, the idea that parallel 
rods can generate a turbulent field 10 - 40 mesh lengths downstream, homo- 
geneous for the purposes of experiments in stagnation regions^ must be con- 
sidered a myth. 

6h Hassler's Averaged Quantitative Field 

In the face of the unsteadiness exhibited in Figs 9 and 10, Hassler 
averaged the signals over one-second intervals (trace lengths over six times 
those in Figs 9 and 10) and labeled this averaged velocity component normal 
to the hot wire u . Fig 13 displays the spanwise variations of u at 
>f = 4* , X'^0.41 cm~2.3S. A vertical cut at x = 2.3^ in Fig 1 tells us 
that u partakes more and more of than of ax-U^ (Equation 11) as y 

increases. Note also that /if a ~0.42.aS for Hiemenz flow and that accord- 
ing to Fig 17 the maximum of is very small and close to the wall at 
K) ■'x/ 1 . 2 or y -w' .09 mm . Thus Hassler's averaging indicates vortical motion 
with A ^ 7+ mm, k~ 0.65 (using 2.4 j'^/a = 1.8 mm), below 

y/^3.6mm= 2^ or r7'-'4.8. While Hassler could document only the peak-to- 
valley spanwise variation of the mean velocity normal to the wire, denoted by 
u*, we can safely infer the mean pattern of motion without direct W^ informa- 
tion because of Colak-Antic 's visual evidence under essentially the same 
conditions in the same laboratory. 

The quantitatively determined mean strength of the vortices , measured 
in terms of u ' , is substantial; see Fig 14. For 'f = 4^ the maximum of the 
ratio u'/u (y=10inm) is approximately 0.2 (0.1 in the U^/U^ , notation of 

Equation (11)). A ’’zero-length" hot wire would indicate still sharper 
minima and maxima in Fig 13 and therefore higher u' . The "instantaneous 
vortices" seen by Colak-Antic must have then been much stronger to yield 
Hassler's average numbers despite their spanwise shifting and dancing. 
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6i The Eactent of t he Hj-e menz and the Gortler-Hammerlin Domains 
and Decay in jx 

The reader who wishes to develop a more concrete feeling for the nature 
of the imbedded Hiemenz-flow approximation will enjoy the exercise in inter- 
preting the mean "normal" u profiles at Y = 4° and in Figs. 13 and 16, 
respectively , in terms of the theory; i.e.. Equations (1) and (2) and the field 
of Fig, 1 when modified for the stagnation region of a circular cylinder by 
Equation (5) . 

For y < 1 mm in Fig. 14, the signal u' should reflect mostly the behavior 
of the theoretically increasing component ax.U^ with downstream distance or Y* , 
but it decreases instead. Presumably this is associated with the curvature of 
the streamlines which now become convex: we are apparently beyond the domain 

II II 

of validity of the Hiemenz and Gortler-Hammerlin flows. The cooling of the hot 
wire by the component may have also contributed to this trend. Note, however, 

that the decay of u' at y = 0.5 mm in Fig. 14 is much less rapid than at y = 2 mm, 
the height monitored in Fig. 15 as a function of y 

Fig. 15 tells us that spanwise variations, u' , decay more or less hand 
in hand with the temporal fluctuations of Fig 10 as the flow accelerates 
around the cylinder. While the choice of the monitoring height may exaggerate 
the decay, it is unlikely that strong strearawise vortical activity would take 
place below this height. The time lines in Colak-Antic' Figs. 11 and 12 
which were generated at the same height of 2 mm at y = 11° , may be displaced 

somewhat inward by the base velocity at this early station, but the main 
upwelling/1 motion in these figures probably occurs near y = 2 mm or even 
farther out. The genuineness of the decay is also born out by mean heat- 
transfer measurements and by the refusal of the boundary layer to become 
turbulent for at Reynolds numbers below on the order of 3.10^ . 

At first glance Figs 11 and 12 of Colak-Antic seem to contradict the decay of 
Fig. 15, but his time-line distortions probably represent primarily the nearly 
frozen integrated effect of past history. It is difficult to judge the local 
time-rate of the distortion which is proportional to the local strength of 
the motion. Colak-Antic agrees with this interpretation. 
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6j The Mysterious Wavelength Doijbling at the Stagnation Line 

The lowest curve depicting variations in the stagnation plane in Fig 
16 seems to show an effect opposite to the tendency to coalescence of vortex 
pairs observed in Fig 12 for larger values. The development of this new 
mysterious effect as the flow approaches the cylinder in the stagnation 
plane is documented in intriguing detail in Fig 16, At y = 10 mm 5 .6 ^ , 
mean hot '-wire measurements (which now correspond purely to - /•’'a^CV^+V^s C-n kz) ) 
sense only the same pure wavelength, manifest at = 4® in Fig 13, but 
somewhat more irregular. At y = 0.5 mm an irregular but definite and 
puzzling double structure is evident — the explanation for the four points at 
f = 0^* marked with a question mark in Fig 21, transcribed from Tani (1974) . 

How can this near-doubling be viewed in consonance with the rest of 
the concepts described here? It could be a manifestation of a true discrete 
eigenvalue selection in a borderline situation where two modes with , say, 

and = 2 would be near equally likely and the randomness of the free- 

stream disturbances would shift the flow between the two patterns in time. 
Since the fields associated with the smaller-scale vortices also reach out 
less in the y direction, it would be possible for the time-averaged 
measurements to sense the larger A^ mode alone at larger distances such as 
y = 10 mm. While not impossible, the situation has no counterpart among 
linear shear-flow instability systems, nor among prototypes of nonlinear 
systems such as discussed in Sections 3 and 4 , 

An alternative is the coexistence at all times of larger eddy formations 
for large y values with smaller formations near the wall. This pattern is 
even harder to imagine and believe. Presumably the mean streamline pattern 
would have to be somehow generalized from that in upper Fig 20 to accommodate 
first a row of free stagnation points at distances apart at a height 
near 7 mm. Closer to the wall a second set of free stagnation points would 
apparently be called for. The two sets would have to mesh in such a way as 
to satisfy the topological consistency of the streamlines and achieve stable 
equilibrium everywhere. Recalling the discussion in Section 3c one tends to 
be pessimistic. 

Any explanation must not only clarify how the mean vortices could 
suddenly split in midstream as the flow approaches the stagnation line but 
also how they could suddenly disappear as the flow leaves the stagnation 
line in the x direction . According to Fig 13, at = 4°, i.e. x = 4.1 mm 
= 2.3 ^ , the smaller vortices are no longer sensed by the probe. 
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In the opinion of the writer, the most probable explanation is that 
the appearance of the second set is induced by the additional boundary con- 
ditions due to the presence of the probe. Stagnation regions and nearly 
separating (i,e,, nearly stagnant) boundary layers are notoriously sensitive 
to probes approaching the wall. Here in the stagnation plane the absence 
(or smallness) of U to relieve the situation compounds the danger of 
localized probe interference. As the probe is shifted to 'f = 4® where its 
momentum defect is no longer convected toward the stagnation line, the 
vortex formations in the stagnation plane cease to be influenced and the 
false signals disappear. In a 1977 lengthy private review of the likelihood 
of the various possibilities, Dr. Hassler also agreed that the last one 
appears as the most probable at this stage of our research. The question 
mark was placed in Fig. 21 to warn of the probability that only the wave- 
lengths reported for 'f= 4^ had best be used for comparison with theoreti- 
cal modeling. 

Another explanation for the measurements of Fig. 16, was pointed out 
by Ivan Beckwith and Dennis Bushnell, of NASA Langley, after reviewing a 
draft of this report. In presence of the vortices spaced A apart, a subset 
of counter rotating vortices with a wave length A/2 may be induced by 
viscous effects and could exist beneath the most unstable mode of vortices. 
Due to the structure of the mean flow (Fig. 1) this subset of vortices would 
move closer to the wall with increasing , so that at ^ = 4° they would 
presumably be found only for y < 0.5 mm. Such vortices have been predicted 
by numerical solutions of the Navier-Stokes equations and observed experi- 
mentally in flows inside cavities beneath a shear layer or near horseshoe 
vortices upstream of an obstacle immersed in a boundary layer. This 
explanation is somewhat related to the second of those listed above. 
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7 OVERVIEW OF ISSUES 


7a Is There a True Instability? 

Perhaps the best point of departure for this Section are the theoretical 
findings of Wilson and Gladwell (1978) . Having rejected the algebraic solutions 
utilized by all others (e.g. in Figs 7, 17, and 18) , they show numerically that 
that Gortler ecjuations then yield only stable solutions. They add, "There is, 
however, considerable experimental evidence that instability of this kind can 
occur," and "It must be concluded that this instability is, as yet, not satis- 
factorily explained" . 

Is it clear that the observed phenomena really constitute an instability 
in the usual fluid-mechanical sense where a linearized system of equations des- 
cribes the basic mechanisms and admits exponential growth in time or space? Do 
we have a well-posed mathematical problem? Wilson and Gladwell tellingly point 
to the severely limited upstream feedback from the body boundary layer and yet 
find it "appropriate to assimie that all distxirbances are proportional to 
exp^ 7‘)exp('i k " • Our Section 4b on Taylor vortices described the conditions 
physically appropriate to exp^T growth with ^ real : disturbances are not 

convected away but reinforce themselves in the closed domain. What physical 
agent could make all the fluid disturbances grow in unison exponentially in time^ 
even far upstream , in otur problem? Experimental evidence illustrated in Fig 8 
suggests that the upstream influence is indeed very limited. Vorticity pertur- 
bations are convected in, interact, and are convected out. Rapid variation of 
the base flow in the direction of convection also distinguishes our flow from 
quasi-parallel flows where exp/BT" (with^ complex) may be rationalized. It was 
the purpose of the studies of the three flows with streamwise vortices in 
Sections 3 and 4 to compare the actual instability characteristics which were 
originally considered analogous. Sections 4a and 4d, in particular, underscore 
the special character of the stagnation flows which make it unlikely that 
6xp/3?- is appropriate. The experiments of Hodson and Nagib in Section 3e, 
however, did not definitely exclude exp^ ? as a possibility. 

7b Implications of Controlled Steady Experiments of Hodson and Nagib 

The same experiments show that steady flows corresponding to the steady 
free stream forcing portrayed in Fig 7 exist. When is small, the forced and 
free response (if any) remains below the resolution limit of the visualization 
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experiments. As is gradually increased, a threshold- like behavior leads to 

the appearance of mature nonlinear vortices with X set by the forcing distur- 
bances in these experiments with limited A range. The onset of the steady 
vortices evidently corresponds to a mild but definite increase in the heat trans- 
fer rate in the stagnation region (Section 6c) . In accordance with Generalized 
Guideline 4 (Section 1) it would be desirable to replicate and improve the 
Nagib-Hodson-Roadman experiments and focus on the stability questions rather than 
heat transfer modules; (see Section 3e for relevant issues) . In particular, the 
widest range of the ratio A / ^ should be tried to uncover any special scale 
receptivity of the stagnation boundary layer as Bippes has done for Gortler vor- 
tices (Appendix la) and Koschmieder for Taylor vortices (Section 4b) , 

7c Forced Motion, Free Response, and Linkage 

Given the nonlinearity of the observed phenomena, can we recognize a 
pattern of known instability behavior in these steady experiments? In Section 4b 
the author proposed a forcing experiment in the non-unique Taylor-cell regime 
using a spanwise arrangement of vortex generators withdrawable through the wall 
of the fixed external cylinder. Withdrawal of the generators after the nonlinear 
cells with admissible A's are established will remove only the steady 
forced contribution to the flow but not the stable cells themselves , Withdrawal 
of the vorticity forcing rods in our stagnation flow removes the whole pheno- 
menon. Forcing and receptivity are central to our problem. 

In the Gortler instability in concave boundary layers and in the Tollmien- 
Schlichting instability, the receptivity in these open systems includes the diffi- 
cult subproblem of penetration of the disturbances into the shear layers . Its 
resolution for rotational disturbances may require a solution of an auxiliary set 
of equations. In the stagnation flow problem the distirrbances are convected 
directly into the boundary layer along the streamlines of Fig 1 and this process 
is govern ed by t he same differential equations as the postulated instability. In 
the steady linearized case, this leads to the Gortler equations with /3 =0. 

If we could resolve the weak-disturbance cases of Hodson-Nagib experimentally, 
how much forced motion and how much free instability response would we observe? 

In the TS and Gortler instabilities above, the experimentally observed modes, 
presumably indistinguishable from normal modes, rise with the disturbance ampli- 
tude. In controlled acoustic excitation of a free mixing layer, the stimulated 
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nonnal-inode vorticity is linearly proportional to the forcing acoustic pressure 
gradient around the point of separation (before nonlinearity sets in) . Does a 
similar linkage exist in our problem in the linear regime (if there is a linear 
instability) ? 

7d A Crucial Limitation of the Hiemenz Plane 

At first thought the direct-inflow aspect of our problem should make a 
theoretical answer to the question far easier than in the other systems (Section 
5b) . The penetration subproblem is not a problem. But the search for an answer 
is thwarted by the limitations of the formulation of the problem in the Hiemenz 
plane alone (Section 5h) , If and 0 as r) —*■ oo , then 

0 by continuity and so does for the algebraically decaying solution, 

which should be linked to vorticity incoming from upstream. Thus treatment of the 
problem without using at least two asymptotically matched regions precludes 
answers to our quantitative search for a forced solution and a linked^ free 
response (if any) , 

The threshold behavior in the Hodson-Nagib experiments indicates that 
interesting answers could be sought through nonlinear formulation of the problem, 
which should then also include other higher order effects such as body curvature 
effect, etc. The initiative of Tani (1974) and lida (1978a, b) in this direction 
is therefore interesting. However, this effort, too, is hampered by confinement 
to the Hiemenz scales . 

7e Continuous Spectrum and Optimum S 

If the algebraically decaying eigensolutions are admitted, a continuous 
spectrum of eigenfunctions for 0 < k 1 exists. Many of their "physical" 
characteristics for the "neutral" or steady solutions, = 0 f are discussed in 
Appendix 4: "Miscellaneous Observations on the Gortler Model". Since, uniquely, 

the eigenvariable is colinear with the motion, the increasing functional values 
as >7 decreases (say, for COx in Fig 17) can be viewed as spatial amplifica- 
tion. Much can be learned about the complex motion by studying the relative 
shapes of the U^, eigenfunctions as dependent 

on the wave number k. It would be interesting to know more about energy 
transfer in these flows, in particular production and dissipation, which are 
now computable . 
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The presence of the continuous spectrum is most likely due to the absence 
of a characteristic length in the Hiemenz base flow (Section 2) . Extension to 
two or more asymptotically matched regions would probably supply an extra condi- 
tion and thereby criteria for an optimum A . At least two sets of researchers 
are working along these lines at this writing . The path is strewn with traps ; 
the matching should be rigorous . A critique of one such attempt by Kestin and 
Wood (1970) can be found in Wilson and Gladwell (1978) . Tani (1974) also pointed 
out that the crucial assimiption of Kestin and Wood to the effect that the vorti- 
city *7 ^ zeros in the whole ^ domain is incorrect (see Appen- 

dix 4) . There has been no rebuttal to the latter criticism, lida (1978a) 
proposes two weak criteria for A selection based on integral properties of the 
eigenfunctions in the Hiemenz domain, 

7f Agents of Instability 

Focusing on the instabilities associated with the vorticity in the bovind- 
ary layer (rather than on the vorticity coming in through the open streamtubes) . 
Wilson and Gladwell (1978) diminish the degree of freedom among the independent 
solutions of the Gortler equation by rejecting the algebraically decaying solu- 
tion as unrepresentative of the possible instability mechanisms. They show 
numerically that for 0 < k <10 the least negative ^ in the exp^ T 
factor is -1,3754 , whatever such global decay implies physically (Section 7a) . 
Having removed the continuous spectrum and believing that the short upstream 
influence is fully covered in the Hiemenz range, Wilson and Gladwell see no 
reason for using matching and inner Hiemenz and outer, y/D , scales, 

advocated above for different purposes. 

They also demonstrate that centrifugal instability is not a factor in 
stagnation- flow instability despite frequent invocation of the analogy in the 
literature. Examination of the terms in the vorticity equation (Section 5a-5c 
and Appendix 2) identifies the stretching of any x - oriented vorticity, as it is 
convected toward the body, as the dominant destabilization in the steady motion. 

It is the Brown University group, Sutera, Kestin, Maeder, Sadeh, Williams, Wood, 
etc. who emphasized that this mechanism (so familiar in turbulence) is at work 
here where others invoked centrifugal instability. However, it was not made 
clear whether they viewed their eigensolutions as forced motion, free normal 
modes, or both. This could make a difference in starting the computations at 
large values — see the paragraph "Computer Differences?" in Appendix 4. 
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7g Experimental Information in Unsteady Environments 

The preceding issues dealt with steady or exponentially growing behavior 
in presence of Fourier-analyzed steady spatial nonhomogeneities in the spanwise 
z direction. Most free-stream environments are also temporally nonhomogeneous , 
in fact more or less turbulent. The Hodson-Nagib experiments with controlled 
onset of Karman instability in the upstream forcing wakes (Sections 6b and 6c) 
provide some insight into sporadic, more violent vortex interactions which are 
possible even for relatively low unsteadiness levels (e.g. Fig 8), The 
careful stereoscopic hydrogen-bubble visualization of Colak-Antic for mild, 
presumably isotropic and nominally homogeneous turbulence add more physical 
insight (Section 6f and Figs 11 and 12). Hassler's unique hot-wire information 
at slightly higher Tu 0.85% (Sections 6d, 6e, 6g, and 6h) brings forth two 
questions : what mechanisms select the dominant temporal quasi -periodicity 

(Pigs 9 and 10) and the prominent average spatial quasi-periodicity 
(Figs 18 and 20) at 240 meshlengths downstream from a nominally homogeneous 
turbulence screen? These most detailed of the available experiments, carried out 
in a reasonably controlled and monitored environment, give indications that the 
spatial vortex structure scales to some extent with the thickness of the boundary 
layer S , Fig 21, The grosser experiments of Brun et al (1966) and Kestin and 
Wood (1970) with less assurance of homogeneity are in general agreement with the 

V V 

results of Hassler and Colak-Antic concerning the one feature of preferred 
average A , Fig 21 . 

At higher intensities of turbulence with distinct mean nonhomogeneity, 

Sadeh and Brauer (1978) demonstrated visually how complex and violent the vortex 
interactions can be. There is little chance to model these types of highly non- 
linear unsteady nonhomogeneous turbulent flows rationally for some time so that 
the conceptual issues perhaps should be addressed in flows with Tu 1%. 

7h Interpretation of Experiments with Unsteady and Turbulent Free Stream 

Do these motions represent forced response, true instability, or a combi- 
nation of both? If we cannot reach an answer to these questions in the case of 
steady disturbances, we shall not be able to do so for this more complex motion. 
The limitations of the confinement in the Hiemenz plane are undoubtedly just as 
operative. We can certainly say from the visualizations that forced motions are 
presen t, probably to an increased degree; see, for instance, the strong forced 

. . '' . V 

vortex xnteraction in Fig 8. Colak-Antic and Hassler agree that in their case 
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the disturbed motion subsides as the vortices are stretched around the cylinder 
with increasing y (Figs 10, 14, and 15) . Such damping is often observed in 
thin laminar bo\mdary layers despite the continued vorticity stretching which, 
however, occurs in the proximity of the wall. Time-line and streakline visuali- 
zation portrays the integrated history of the fluid elements , not their local 
rates of deformation, so that Figures like 12 and the stills of Sadeh and Brauer 
(197fci) often give false impressions of excessive activity. Even at relatively 
high body Reynolds numbers the agitated motion does not cause turbulent bursts 
intrinsic to the scales of the body boundary layer. Thus the enhanced heat 
transfer subsides around the cylinder \intil a genuine transition in the buffeted 
laminar boundary layer of the body takes place at high Reynolds numbers . The simplest 
telltale sign of transition is the onset of high frequency oscillations in the botindary 
layer on the scale of 0.1 ^ to 0.0015 and less; (see also Appendix 3). 

Does the evidence in Fig 21 of vortex scales on the order of 2 - 4 5^ 
prove that there exists genuine laminar instability in the stagnation region in 
the sense of Section 7a, alongside the forcing motion? It suggests the possi- 
bility of unstable response but by no means demonstrates its presence. The forced 
response, after all, is also a solution of the motion in presence of the boundary- 
layer field and must reflect some of its characteristics such as the thickness. 
(Similarly the frequencies in decaying Wygnanski -Champagne puffs in a pipe at 
Reynolds numbers below 1800 will scale with the pipe shear layer without consti- 
tuting true instabilities of the shear layer.) 

7i Whither Further Progress? 

This brings us back to the question we posed in Section 7a in connection 
with the apparent impasse of the Wilson-Gladwell findings: is there a true 

instability in the steady and unsteady cases? The question is unlikely to be 
answered experimentally because the instruments read the whole motion, and it is 
impossible to decompose the signals without knowing much about the characteristics 
of the ingredients. Theoretically the most promising approach seems to be the search 
for the linkage between the forced and free responses (if any) for the steady 
environment, which means breaking out of the Hiemenz scales constraint 
(Sections 7c - 7d) . 

For the unsteady problem it might be worthwhile to explore changes in the 
structure of the linear problem (forced or free) which would ensue if one replaced 
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exg^l' { real) with ©3^ ^i3T7‘) ( iT real or complex) in Gortler's 
formulation, Equation (11), Kayalar (1969 ) set up some general equations along 
this line but studied only the case with k=0 , i.e, the bulk oscillation 
without vortex formation. Could there perhaps be any new ways of transferring 
energy to the distiorbed flow due to special phase relationships? 

Since the Hassler hot-wire information is currently unique, it should 
serve as a springboard for still more incisive and more controlled experiments, 
in accordance with the Generalized Guideline No, 4 of Section 1. If two-probe 
space-time correlations could be accommodated for this interference-prone 
geometry, the experiments should prove most rewarding. 

Without improved information on detailed motion and on coherent motion we 
can say little about correspondence of unsteady reality with theoretical models . 
For instance, to use any comparison of optimal A values^ drawn from steady 

linear theory in the Hiemenz plane, with the measurements in Fig 21 in unsteady, 
highly nonlinear flows as sufficient evidence for a model is an exercise in 
lanwarranted extreme faith. Analyzing the present experimental information makes 
one feel very humble and conscious of the need for rigorous experimental proce- 
dures and rigorous inferences . 

7j Applicability of the Information 

Although the original motivation for scrutinizing the stagnation flows 
came from applications, the emphasis in this report has been on the basic infor- 
mation. Reliable, rational semi-empirical prediction methods for a host of 
applications cannot be expected without deeper appreciation of the physics and 
mathematics of the problems . 

The most immediate applications deal with heat and mass transfer in stag- 
nation regions of bluff bodies. When the leading edge is quasi-cylindrical and 
normal to the flow, the present account of the experiments and theory provides 
us with some conceptual understanding of the underlying mechanisms . It should 
appeal to the intuition of the engineer. Both the vortex roll-ups associated 
with spatial nonhomogeneity in the environment and the unsteady vortex inter- 
actions in a homogeneous turbulent environment contribute to an increase in 
scalar property transfer (Sections 6a - 6c) , without the body boundary layer 
becoming turbulent. Unless the body boundary layer undergoes transition to its 
own intrinsic turbulence, the heat transfer enhancement subsides with increasing 
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<f . The motions associated with the enhanced heat transfer evidently bear a 
resemblance to the general pattern of motion envisaged by Gortler in 1955 » The 
vorticity -stretching mechanism is part of the dynamics, as has also been recog- 
by Bearman (1972) cind Hunt (1973) , 

The subject of transfer mechanisms and boundary layer instabilities at 
swept-back cylinders auid axisymmetric bodies is a separate matter. There seem 
to be some similarities cind some important differences so that drawing analogies 
is not recommended without further studies. 
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APPENDIX 1 


RECEPTIVITY, INITIAL VALUE PROBLEMS, TEMPORAL AND SPATIAL NORMAL-MODE 

APPROACH AND NONLINEARITY 

a. The Problem of Receptivity to Disturbances and Som e Illustrations 

II ' 

In a detailed study of the growth of Gortler streamwise vortices in 

boundary layers on a concave surface in water Bippes (1972) deliberately 

stimulated disturbances of different spanwise wavelengths X by local 

heating at the wall. When the imposed X was substantially larger than 

the eigenwave length X^ of maximal amplification the initially doninant 

flew response was in vortical motions with the forcing periodicity X, 

followed by slow spatial development of ultimately dominant motions near 

X^. When Knapp and Roache (1968) irradiated the laminar boundary layer 

on an axisymmetric ogive-cylinder body with an acoustic frequency f^ in 

a relatively broad band near the maximally excitable Tollmien-Schlichting 

frequency f^, the response apparently developed into that of free TS eigen- 

waves at the forcing frequency f^. When loud sound at frequency f^ near 

f^/2 was beamed at the boundary layer, the response appeared to be a TS 

wave of frequency 2f^, again near f^. When Loehrke (1970) gave a sudden 

pulse to a spanwise stretched thin wire (0.038mm in diameter) , vibrating 

in the flew direction within a 4mm thick Blasius boundary layer, a slowly 

broadening wave packet of TS waves with frequencies around f grew and 

m 

traveled downstream. 

These are but three examples of selective responses of unstable 
shear layers to identified forcing disturbances. Complicated responses 
with wavenumber spectra in the region of most-amplified eigenfuctions 
are knewn to evolve in presence of lew-intensity turbulence in the free 

If 

stream for both the TS and Gortler instabilities. The intensity of the 
response grews with the intensity of the free-stream disturbances and 
leads to earlier transition to turbulence in the shear layer. In open 
flow systems: duct flows, free shear layers, and boundary layers (in 

wind tunnels or in free flight) , the upstream regions generally contain 
nonvanishing acoustic, vorticity or entropy fluctuations which penetrate 
the shear layers and excite their responses. These responses evidently 
include free eigenfunction motions (normal modes) which dominate in vin- 
stable shear layers. The processes of penetration, excitation of forced 
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and free responses ("internalization” of the disturbances) constitute the 
receptivity problem as identified by Morkovin (1969) cind discussed by 
Reshotko (1976) . Rigorous experimental elucidation of receptivity and 
its convincing theoretical modeling represent an important step in our 
further understanding and estimation of transition to turbulence, see 
Morkovin (1978) , 

b. Temporal Initial Value Problems 

Receptivity is sometimes linked to initial value problems of linear 
instability. For such problem formulation it is assumed that the dis- 
turbances have already penetrated the shear layer and that they are knom 
throughout the region at some initial time t^. The evolution of the forced 
and free responses is sought for generally in the easier temporal 

formulation. The motion associated with any normal mode k then grows 
(or decays) uniformly throughout space in proportion to exp(3j^t) . This 
approach has the virtue of potentially identifying the forced response 
which presumably has little to do with the build up of disturbance energy 
to nonlinear levels. (However, this motion may possibly play a role in the 
initiation of subsequent secondciry instabilities which seem to be necessary, 
but not necessarily sufficient, for onset of turbulence) . Quantitative 
linking of the disturbance at t^ with initial amplitudes at of the 
amplified normal modes of the free response would be conceptually helpful 
to the understanding of one element of receptivity. The equally important 
element of receptivity, associated with the transformations of the distur- 
bances as they penetrate the shear layer and adjust to it, is absent in 
the initial value problems. 

c. Spatial Initial Value Problems 

We must keep in mind that the analytic problem is exceedingly diffi- 
cult and that without substantial idealization it becomes intractable. 

But whatever the useful conceptual lessons of initial value problems 
in time are, it should also be recognized that the approach does not seem 
to be adaptable to a nxamber of important situations. For instance it 
does not fit the cited cases of Loehrke , or Knapp and Roache, and of 

II 

Bippes, nor Bippes* case of excitation of Gortler instability by free- 
stream turbulence. In the cases of Loehrke and Bippes, the disturbances 
were introduced artificially directly into the boundary layer at fixed 
streamwise locations, and the forced cind free responses evolved spatially 
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rather than in time. Apparently, the single-frequency, constant-amplitude 
acoustic excitation of Knapp and Roache also developed according to the 
spatial instability formulation, but involved furthermore the penetration 
and internalization problems. For these cases an initial-value problem 
in the streamwise direction x would be more appropriate although it is 
mathematically less well posed. (The handling of the downstream boundary 
conditions is not rigorous and remains hoc at this writing) , Numer- 
ical experimentation on conputers with either linearized or full Navier- 
Stokes equations, e.g. Fasel (1976) , Fasel etal (1977) and Orszag (1979) , 
suggests that additional insight into receptivity could be gained with this 
approach for two-dimensional and simple three-dimensional disturbances, 
given enough funds and computer capacity. The word experimentation is 
used on purpose to emphasize that one loses the generality of analysis 
and that there are also problems of resolution and convergence. On the 
other hand nonlinearity can be handled. At this writing, however, no pub- 
lished contributions to the receptivity problem have appeared based on the 
spatial initial value problem. 

d. Issues in Spatial and Temporal Growths in Normal- Modes Techniques 

In principle, the temporal or spatial initial value problems can provide 
a quantitative sequential link between the internalized initial disturbances 
and the free normal modes which ultimately dominate the unstable motion. 

However, most linear instability analyses merely assxime exponential temporal 
or spatial growths of normal modes of specific functional form and focus on 
the analytic-nxamerical problem of determining the eigenvalues and some typical 
eigenfunctions. Again the assumed temporal or spatial formulation may not fit 
the physics of the problem and may lead to serious difficulties in interpreta- 

II 

tion. Thus, in our first example of receptivity, the streamwise Gortler vortex 
instability was originally assumed to grow as exp3t while Bippes (1972) and 
others find only growth in x. For Tollmien-Schlichting propagating instabilities 
Caster's transformations establish a direct link between spatial and temporal 
growth of the normal modes which is valid for relatively slow amplification rates 
in non-decelerating quasi-parallel boundary layers. Caster (1962) . 

For a system to be capable of both exponential growth in t and x, 
implicit in the Caster transformation, the linearized stability equations 
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must have coefficients independent of both t and x, i.e. they should 
correspond to a strictly steady and parallel base flow of infinite extent. 

For slowly growing boundary layers the slow x variation leads to second - 
order corrections which do not vitiate the physical interpretation. 

For non-parallel base flows with rapidly converging or diverging stream- 
lines, such as our flow toward the stagnation line, the exponential 
disturbance growth in the streamwise direction is precluded by the variable 
coefficients in the governing equations. Since the coefficients remain 
constant in t, the expBt factor for the solutions is available - and used 
by default. Its meaning deserves careful scrutiny. 

For quasi-parallel boundary - layer instabilities of the TS and 

II 

Gortler type, the eigenvalue problems develop in the cross-stream y 
direction. Free-stream vorticity and entropy disturbances travel along 
mean streamlines of the base flow and influence the fluctuations within 
the shear layer either by vortical Biot-Savart induction across streamlines 
or by slow penetration along with the external streamlines because of the 
departure from strict parallelism. (Acoustic disturbances, of course, 
can always propagate across streamlines) . 

e . Some Special Features of the Eigenvalue Problem for Stagnation Flows 

Here again the contrast with the stagnation flow near the plane of 
symmetry x - 0 in Fig. 1 is striking. All vorticity and entropy dis- 
turbances are directly convected into the region of presumed instability. 
Furthermore, unlike in the boundary-layer cases, they are subjected to a 
known straining field (due to the base flow) during this convection. This 
means that forced response of the flow to upstream disturbances should 
be easier to calculate. However, the singular fact that the independent 
variable y or ^ in the eigenvalue problem is here in the direction of the 
stream and thus ties together the whole upstream history of the flew 
suggests that the physics of the eigenrelations is basically different. 
Geonetrically , this is associated with the streamline bifurcation and 
the 90% turn at the stagnation line, not present for other instabilities. 

Sensing the peculiarity of the variable^ Wilson and Gladwell (1978) 
argue on p. 517 that "the instability, if any, is the centrifugal insta- 
bility" which can take place only as the concave streamlines reach the 
vortical region of the base flow, i.e. the viscous layer. They believe 
that the upstream influence (working against the inflow) of a mechanism 


63 


confined to the viscous layer should be dying out exponentially and 
thereby they change the character of the eigenvalue problem. 

Wilson and Gladwell consider the possible role of free-stream 
disturbances and the effect of the upstream straining only in their 
closing remarks and negatively so. A major part of their argument is 
based on direct misreading of the German text of Hassler (1971) and 

^ V 

Colak-Antic (1971) . Their experiments did not have "exceedingly small 
free-stream vorticity" because they towed a turbulence-generating grid 
ahead of their cylinder in the otherwise quiet water tank- The reported 
measured intensity is 0.85%, i.e. substantial enough to be dominant in 
other instabilities . The other part of the Wilson-Gladwell reasoning is 
taken up in the final discussion of the phenomenon. 

f . Reflections on Receptivity to Turbulence Leading beyond Initial 
Value Problems 

Since time and again interpretation of instability experiments hinges 
on the role of free-stream turbulence, it is worth-while to reflect on the 
nature of such turbulence and on how the inherent variability of its local 
intensity and scale may interact with potentially exponential responses in 
shear layers. It appears likely to this observer that receptivity to free- 
stream turbulence involves a problem more complex than can be handled by 
either of the initial value problem formulations. Already the first 
experiments on TS waves, those of Sch\ibauer and Skramstad (1947), disclosed 
that free-stream turbulence elicits modulated wave-packet response in 
boundary layers . Free-stream turbulence in uniform flows is a decaying 
field, spotty in space and time (and only nominally isotropic and homo- 
geneous when time or ensemble averaged) . Especially in the latter stages 
of decay, there are sporadic occurrences of higher-amplitude events and 
of patches with different scales, occasionally referred to as inner intermittency . 
A boundary-layer response in the foimi of finite TS wave packets is therefore 
consistent with such spotty, time-modulated, three dimensional, forcing 
disturbances. Within a given wave packet there is energy transfer in space 
and time. Because of this energy sharing the primary packets induced in the 
upstream region of a growing boundary layer by a spotty sporadic vorticity 
formation as it is convected downstream near and within the boundary layer 
will probably have to be characterized by a finite spacetime domain, Ax, Az, At. 
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We note that information on temporal averages of free-stream 
turbulence or even cai its three-dimensional spectra (if obtainable) 
would be of little help for the above picture of response to weak 
turbulence. Rather, for fundamental studies of receptivity at least 
two-probe, space-time correlations of the exciting and response signals 
are indicated . 

g. Apparent Consequences of Substantial Amplitude Variations in 

X and t: Gaster Experiments 

The exploration by Gaster and Grant (1975) of simpler wave packets 
generated by sudden puffs through a small hole in a plate supporting a 
growing Blasius boundary layer and its theoretical modeling, Gaster (1975) , 
taught us much about the linearized development of three-dimensional TS 
wave packets. It appears now, Gaster (1978) (subject to Gen. Reshotko 4), 
that the finiteness of the wave packet in A z. At, brings about 
"dramatically" different nonlinear development and an earlier secondary 
instability (also called breakdown of TS waves) . This comes from compar- 
isons of the behavior of the wave excited by the sudden puffs and those 
excited by "continuous puffing" through the same hole at the central 
frequency f^ of the pulsed wave packets. Strictly speaking this is not 
a receptivity effect, because the linear development (which is present 
because the disturbance is linearizably small) is modeled satisfactorily. 
However the effect warns us that there may be dangers of misinterpretation 
when we use the otherwise attractive procedure of studying responses to 
single-frequency normal modes. 

h. Subcritical Instability, Threshold Phenomena and Receptivity ; 

The Duct-Flow Case 

When the environment disturbances are la,rger and involve some non- 
linear processes the problem of receptivity gets even more complicated. 

It is likely that the cited frequency response of the boundary layer 
of Knapp and Roache to a forcing acoustic disturbance f^ resulted from 
nonlinearity in the sound. The harmonic 2f.^ component may have been 
hidden in the background noise and yet brought out by the exponential 
amplification of the TS mechanism. Analysis of weakly nonlinear, slightly 
stable systems, e.g. Stuart (1971) p. 359, shows that "nonlinear effects 
may permit the existence of a threshold amplitude" above which instability 
can grow even though the system parameter is below the critical value for 
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linear Cinfinitesimal) instability. 

Hard evidence of such subcritical instabilities may elude our measuring 
instruments especially if it develops temporally in a truly parallel flow. 
For instance it has been stated in niamerous papers that TS waves in two- 
dimensional Poseuille flow cannot be observed because the growth of the 
instability would be so rapid that only the final turbulent state of 
equilibrium would be ascertainable. And indeed for more than half a 
century no investigator was successful in preventing transition to turbu- 
lence above 66 to 90 percent of the infinitesimal critical Reynolds number, 

- until Nishioka, lida, and Ichikawa (1975) brought down their free-stream 
turbulence to approximately 0.01% and maintained laminar flow beyond 
Re .see Morkovin (1978, pp 7-9) for discussion of the issues . By intro- 
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ducing controlled, nearly two-dimensional disturbances through a vibrating 
ribbon^ Nishioka etal also set the stage for a spatially growing initial 
value problem and demonstrated the spatial subcritical threshold behavior 
of their TS waves. It is instructive with respect to the character and 
role of disturbances that Karnitz, Potter and Smith (1974) were able to 
reach only 0.87% of with free-stream disturbances of approximately 

0.3% and reported privately that any loud voices or even finger snapping 
near the Poiseuille duct would switch the flow from laminar to turbulent 
near that Reynolds number. Theory has established that this flow exhibits 
subcritical bifurcation so that we are forewarned in our interpretations 
of the influences of disturbances. But what about interpretation of 
behavior in presence of disturbances in systems that have not been adjudged 
subcritically unstable? Or by observers who are not familiar with the 
peculiarities of nonlinear instability? 

i. Diagnostics, Thresholds, and Visualization Evidence 

The awareness of the possiblility of threshold phenomena is an essen- 
tial ingredient of diagnostics of instabilities. The willingness and 
capability of devising deliberate perturbations (the general "spoiler" 
technique, Morkovin 1969) is another mainstay of instability diagnostics. 
From such explorations by Hodson and Nagib (1975) it appears that the 
stagnation flow does show threshold behavior with respect to steady low- 
Reynolds number disturbances, see Section 3d and 3e. 

An important question with respect to diagnostics relying on visuali- 
zation techniques is their capability of resolving disturbances in their 
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linearizable stages. General information from hot-wire anemometry, 
primarily for TS type instabilities, suggests that velocity fluctuations 
below approximately 1% of local maximum velocity exhibit behavior 
consistent with linearized equations, e.g. apparent local exponential 
growth, cibsence of higher harmonics (within resolution of instrument) , 
independent superpos ability of disturbances with different frequencies, 
no threshold with respect to subsequent stabilization etc . With care 
and know-how, monochromatic fluctuation amplitude of the streamwise 
velocity component can be measured quite reliably to 0.01% of the local 
maximum velocity so that the 1% estimate cannot be grossly in error and 
certainly is not on the order of 3% . 

Comparison of various visual indications of spatially growing 
narrow-band TS type waves with hot-wire traverses of the same controlled 
fields lead one to conclude that visualization almost surely does not 
resolve the fluctuations in their linearizable stage,* The conclusion 
is reinforced by studying Bippes' 1972 stereographic measurements of 
steady Gortler vorticies growing in x. The implication of these 
considerations is that when using visualization diagnostics we should 
be conscious that we are probably seeing only the nonlinear regimes^ and 
that a threshold phenomenon may have taken place. Generally, then, we 
should distinguish between linear receptivity following a threshold 
jump when making our diagnoses. 


* Actually the issue is sufficiently important in diagnostics to 

warrant quantitative research for the several visualization techniques 
smoke, hydrogen bubbles, thin oil films which indicate stress lines at 
the wall, evaporation and sublimation techniques. It is possible that 
steady streamwise vortices may be disclosed in their linear stages by 
the cumulative effect of evaporation and sublimation over sufficiently 
long time periods. 
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APPENDIX 2 


ON yORTICITY AND PRESSURE RELATIONS 

For ’’incompressible'' flows with, constant viscosity but variable ( for 
instance stratified) density^ application of the skew curl operation to the 
Navier-Stokes momentum equation, V :x M , leads to an equivalent of a local 
angular momentum equation, expressed in terms of vorticity ^^u. = w 


D(u)/p) 

Dt 


V „2 w 
_ V (D = -— 

p — p 



(A2-1) 


This skew operation is said to decouple vorticity from its companion, the 
pressure, which in the momentiam equation appears as ^p , i.e. operated upon 
by a symmetric operation. 

When we use the symmetric divergence operator on the Navier-Stokes 
momentum equation, > we banish skew entities by virtue of V - u) = 0 

and arrive at the pressure equation, companion of the vorticity equation (A2-1) 


(now with density constant) 
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For compressible flows both the vorticity and the pressure equations get much 

more complicated, but only Equation (A2-2) changes its nature by becoming hyper- 
bolic: the operator on the left hand side of the generalized equation (A2-2) 

becomes the wave operator (1/a ^) , (e.g. Phillips, (I960)), 

2 

where a is the speed of sound squared. It is then possible to use the 
equation as a vehicle for analysis of aerodynamical ly generated sound, with the 
non-homogeneous right- hand -side terms of the wave equation interpreted as 
effective generators, attenuators, and redistributors of sound, generally 
called "sources". As a fluid becomes "more incompressible", its speed of sound, 
a , increases indefinitely and the hyperbolic wave operator reduces to Laplace's 
operator in Equation (A2-2) . Thus for idealized incompressible fluids any 
pressure changes, however caused, act "instantaneously" at any distance and 
come back "instantaneously" from any "reflectors" in the field. This mutual 
instantaneous interaction within the whole field conveys one of the essential 
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properties of the resulting elliptic differential equation. Since it at times 
camouflages causes and effects, it is often helpful to think of slightly com- 
pressible fluids with, large but finite speed of sound. The instantaneous 
action concept is also present in the Biot-Savart law of induction of velocity 
by a "source-like" vortex. Such, action must in fact be transmitted "hand in 
hand with the companion agent, the pressure", and often it is useful to think 
of it as propagating at the speed of sound, e.g. for linearized supersonic 
flows . 

The concept of sources comes from the fruitful study of linear differ- 
ential equations with their right hand side given by the Diract delta or source 

function, such as 

^ ^ - V^P = «(x- Xq . t- t^) , (A2-3) 

a dt 

II - K V^T = «(x- Xq , t- tg) . (A2-4) 


The latter equation is the familiar parabolic heat equation for non-deforming 


(stagnant) fluids or solids where K stands for Kelvin's diffusivity or 
Maxwell's thermometric conductivity, e.g. Carslaw and Jaeger (1947), Chapter I, 
The inquisitive reader will also find a lucid presentation of the power of the 
source concept in Chapters X and XIII of Carslaw and Jaeger. For non-stagnant 
fluids, the local derivative operator 2)/^t is replaced by the material deriv- 
ative operator D/Dt so that the left-hand-side operator in the heat equation 
(A2-4) becomes identical with the operator in the vorticity equation (A2-1) 
if K is replaced by V , the kinematic viscosity. 

It is instructive to linearize equation (A2-1) around a two-dimen- 
sional base flow Ug(x,y) , (x,y) , g ^(x,y) for the x-perturbation component 

of vorticity LJ and put it in the above form: 

- Vgif « _ V ^ ^ (^2-5) 

The intrinsically three-dimensional terms on the right-hand side thus act as 

effective sources of uj . The first term is the U) stretching source due to 

X X 

base- flow straining while the other two terms are due to rotation of and 

onto the cj direction. These two terms, however, are automatically 
BZ X 
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zero for the outer Hiemenz flow where /^y ^ terms 

o HZ 

cancel each other exactly in the GortlerrH'inraerlin flow for which and w 
remain independent of :x . In these cases the vorticity stretching indeed 
dominates the behavior of CJ , 
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APPENDIX 3 


ON TURBUIiENCE CONCEPTS USEFUL 
IN DIAGNOSTICS OF LAMINAR-TURBULENT TRANSITION 

TURBULENCE implies spatially and temporally irregular nonlinear three- 
dimensional vortical motion resulting in transport of momentum (heat or mass ) 
much larger than that due to molecular processes in laminar motion at the same 
mean Reynolds number (Rayleigh number or Taylor nxmiber, etc.) . For point 
measurements turbulent character implies continuous three-dimensional power 
spectra decaying monotonically past some characteristic wave number (frequency) . 
Turbulence may be present locally within large-scale organized motions (nearly 
regular, nearly periodic) . 

TURBULENT- LAMINAR INTERFACES occur at sufficiently high Reynolds nimabers 
as sharp gradients in flow characteristics between irrotational (exterior) flows 
and turbulent shear flows (jets, boundary layers, wakes,, plumes, etc.); also 
between laminar shear flows with regular oriented vorticity and neighboring 
shear regions which had ;indergone transition to irregular three-dimensional 
vortical motion (boundary layers, pipe flows, Couette flows with counter- 
rotating cylinders, etc.). 

TURBULENT SPOTS are formed through isolated localized fine-scale, 
higher-order-instability breakdowns of laminar shear flows. Through the 
ensuing turbulent-laminar interfaces turbulence spreads into adjacent laminar 
neighborhoods and is simultaneously convected with a "mean speed". In boundary 
layers "TRANSITION" begins at the upstream locations of the sporadically forming 
turbulent spots and ends where the spots have grown together leaving no laminar 
lacunas . 

INTERMITTENCY refers to repeated crossings of laminar- turbulent inter - 
faces at a fixed observer point. It is usually characterized by T r the 
fractional time (between 0 and 1) of turbulent behavior at the point. 

In TURBULENT BOUNDARY LAYERS identifiable violent motions called " bursts " 
erupt sporadically near the wall. A disproportionally large portion of turbu- 
lent energy and Reynolds stress is generated in the burst activity , Bursts 
do not take place in boundary layers below some still poorly defined Reynolds 
n\amber Re^ min on pressure gradient, three-dimensionality, etc. 
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Interfaces associated with turbulent boundary layers or spots at low 
Reynolds numbers (jiist past ^6 min ^ diffuse, and it is difficult to 

determine operationally whether a boundary layer is a laminar, randomly 
buffeted boundary layer or a self T^ustaining turbulent boundaiy layer. 

BUFFETED BOUNDARY LAYERS , When upstream turbulence impinges upon a 
boundary layer below ^6 min ' resulting motion has all the specified 

earmarks of turbulence, except that the scales of the motion exceed the 
boundary layer thickness , The mean gradients near the wall are then less 
steep than if the boundary layer itself was intrinsically turbulent. The 
transport of momentum (heat or mass) in such laminar buffeted boundary layers 
is therefore enhanced, i.e. more than if the bo\mdary layer were not buffeted 
and less than if it were turbulent on its own scales . As the buffeted boundary 
layer grows, the transport of momentum (heat or mass) will increase as its own 
transition takes place. Laminar buffeted boundary layers occur frequently at 
leading edges of bodies and are of technological importance, e.g. in turbo- 
machinery. 

One can speculate that in high-speed Couette flow periodically oscilla- 
ting Taylor cells with periodic instabilities in wall boundary layers or in 
shear layers between cells can also be buffeted by neighboring turbulent 
regions . Since turbulence sets in at discrete locations , coexistence of 
laminar and turbulent regions should be expected. The temporal spectrum in 
such a doubly periodic cell near another partly turbulent cell would consist 
of two discrete lines plus a continuous spectrimi as is observed. 
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APPENDIX 4 


MISCELLANEOUS OBSERVATIONS ON THE GORTLER MODEL 

Geometrical configurations often possess preferred vortical formations. 
Gortler’s "spontaneous" pattern. Equations (11), or its forced counterpart 
(say, generated by weak merged wakes from the dotted grid in lower Figure 3) 
certainly describe a family behavior which is present in many physical 
situations but not necessarily all. Kayalar (1969) proposed a different set 
of streamwise vortical patterns, without theoretical or experimental backing. 
Inger (1975) discussed another "forced disturbance model" of streamwise 
vortices, which may be too simplified, e.g. V 2 (ri) = 0 , and saw little 
evidence for their existence. In this section, agreements on Gortler forma- 
tions and the corresponding physical fields are reviewed in detail. 

The Detailed Fields 


The general flow pattern, though to second order in the small quantities/ 
is represented in Figure 7 and upper Figure 20 according to Kestin and 
Wood (1970) . The D variation of the perturbation velocity components 
and the pressure for k = 2/3 in Figure 17 was borrowed from Inger (1974) 
where additional prof iles for k = 0.2, 0.5 and 1 , as well as for cases of 
very cold and very hot cylinder walls can be found. The lower part of 
Figure 20 for x = constant was constructed to show schematically the relation 
between the full flow streamlines in the upper part of Figure 20 and the 
linearized , W-j^ , u)^ profiles of Figure 17. The circles with + and 

- signs locate the local maxima and minima of the linear pressure 
perturbation with the ones near the wall being typically 4-6 times stronger 
than those near the edge of the boundary layer. The arrows show the relative 
magnitude of the physical perturbation velocity (which is opposite 

to V ) according to Equation (11) . 

Nonlinearity 

The W-j^ field and the rotational patterns of lower Figure 20 are directly 

comparable with the upper figure, except for the second order effects present 
in the latter. However, for the full vertical-velocity and pressure picture 
we must superpose their counterparts on the base Hiemenz flow. The amplitudes 
of the linear eigensolutions in Figure 17 are arbitrary as long as they are 
"small". One could normally dare to assume ratios of maxima ^ the boundary 
layer of * ®tc. to be on the order of 0.2 for the linear 

regime and up to 0.4 for second order analysis and still retain good pattern 
resemblance. However at the value of n corresponding to the free stagnation 
point in the Kestin-Wood upper Figure 20, where ^ ~ 0 * ^ 

therefore V is locally a rather massive "disturbance". That such .large 
steady disturbances commonly occur in the real world has been shown by Hodson 
and Nagib (1975), e.g. upper Figure 5. 
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The Pressure Fields 


The flow away from the pressure minima at kC * 0 and 2 tt in lower 
Figure 20 could be puzzling until one makes some estimates of the Hiemenz 
pressure rise from the + circle to the — circle that must be superposed. 

2 2 

According to Equation (6) this Ap~ Cat +) = (+)]^ . Because 

it takes a V perturbation of |'fQC+)| ~5 to cause stagnation at the + 

circle, and because the largely viscous pressure perturbation is likely 

to remain nearly proportional to the first power of the velocity perturbation 

this quadratic inertial pressure rise will surely make the full static 

pressure at the - circle larger than that at the + circle. The pressure 
at the + circle near the wall at z= X/2 will then be the highest 
achieved stagnation pressure. The relative magnitudes of the stagnation 
pressures at the three stagnation points , S and in upper Figure 20 

should then be in the same descending order as at , S and in lower 

Figure 3 discussed in connection with the threshold vortex appearance in 
Section 3e. In the linearized viscous solution (lower Figure 20), the local 
pressure minima seem to be induced by vortex associated entrainment (viscous 
pumping) away from the regions in question. 

Does Vorticity Have a Zero? 


According to lower Figure 20, the standing vortices have cores with solid 
body rotation as witnessed by the linear distribution of velocities from 
zero at their centers. But there are no pressure minima at these centers, 
as there would be in vortices resulting from roll up of thin vorticity sheets 
It is a topological necessity for these vortices to induce velocities 

as shown, velocities which generate their own thin vorticity boundary layers 
at the wall of a sign necessarily opposite to that of the standing vortices. 
Thus the upper figure, borrowed from the 1970 paper of Kestin and Wood, 
demonstrates the impossibility of a key statement of theirs in the same 
paper, to the effect that does not change signs between q = 0 and 

f) “ (we shall return to the implications of the statement in Section 7 ) . 
All linear vorticity computations (including those of Kestin's coworkers 
Sutera (1965), Figure 4, and Williams (1968), Figure 11) disclose the sign 
reversal as in Figure 17 here, but the topological argument transcends 
linearity. 

The Smallness and Limited "Reach" of U Fields 


Thus with minor exceptions, there is a general accord on the consistency 
of the Gortler vorticity and velocity patterns and their occurrence in steady 
real flows, Hodson and Nagib (1975). All the computed perturbation profiles: 
Sutera (1965) , Williams (1968), Inger (1974) and Tani (1974) also agree on 
the shapes of the profiles, their sign changes and on the smallness of the 
amplitude and the limited q extent of the perturbation (Figure 17). 

Sutera and Inger observe that acts as if it were decoupled from 
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and > Since velocities in the y and z directions and the pressure 

in Gortler’s model. Equations (11) are Independent of x , the coupling of 
the equation with the others is solely through the convective derivative 

SUp/By ■' » which vanishes outside of the Hiemenz boundary layer 

near n “ 4 . Thus the equation is indeed independent of the rest, 

outside of the unperturbed boundary layer. Hammerlin (1955) used that fact 
to start his asymptotic analysis. The single independent nondivergent 
asymptotic solution for the second order differential equation in has 

2 *. 2 

the leading term Af(k)exp^ — — 0^648) j ^ 0.648) . The quadratic 

exponential decay signals the quenching of the coupling source. (A is one 
of the free coefficients of the linearly independent solutions discussed in 
Section 5 in connection with Hammerlin' s smooth matching of the regular inner 
solutions and the outer asymptotic solutions.) The extensive outer reach 
of and variations relative to that of underscores the inviscid 

nature of the basic amplifying processes (Equation (9), Section 5 , or 
— (0^3V/8y)i the component is stimulated only through the gradient 

BU^/By which is generated within the viscous boundary layer. 

Contrast with Gortler's Vortices on Concave Walls 


In 1955, before computers could numerically experiment with complex 
equations and grind out results on which the present discussion is partly 
based, Gortler conjectured about the analogy between the vortices in the 
stagnation region and the "other" vortices bearing his name, arising in 
boundary layers on concave walls. The latter are driven by centrifugal 
instability, just like the famous Taylor vortices exterior to a rotating 
cylinder, whereas streamline divergence and adverse pressure gradients are 
essential here. It is unfortunate that this early conjecture is being 
perpetuated in widely used texts and reference books. The structure of the 
"other" Gortler vortices is vastly different . In particular the 

component, instead of being smaller than and by factors on the order 

20 (for reasons already discussed) is orders of magnitude larger as pointed 
out by A.M.O. Smith (1955) and Betchov and Criminale (1967). In one case 
Smith computed the ratios and /Wlmax and 

667, respectively! 

The Spanwlse Flow Field 

The dimensional steady-state velocity in the x direction, 

(h)+ axU^ (H) cos kC , also have x as coefficient, making the 

perturbation term still smaller in the region of validity of the approximation 
The normal component of velocity, ($q+Vj^cos k^) , grows only as 

from the wall because V(0) = V'(0) = 0 . The spanwlse component is by far 
the largest in the vicinity of the wall at x = 0 , In other words, the 
solutions suggest that considerable velocity "can be generated along the 
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stagnation line toward S " despite the strearawlse "leak" promoted by 
acceleration a^j (see the end of Section 3c ) but that this motion does not 

build up a larger stagnation pressure pg because of the viscous pumping 

effect near S . The streamline pattern near the stagnation line should 
therefore look as depicted in Figure 6 in inger (1977) , 

The Temperature Field and its Limited "Reach" 

The slow streamwise removal of matter from the stagnation region should 
also make the steady-flow heat transfer at the wall more limited than 
intuitive ideas comparing molecular "mixing lengths" of the undisturbed 
Hiemenz flow with the "vortex-size mixing lengths" of the perturbed flow 
would at first suggest. Insofar as the fluid chums without substantial 
streamwise motion the heated elements return to the wall decreasing the 
original AT between the fresh fluid and the wall and therefore the 
capacity for subsequent heat flow. In the steady nonlinear regime there is 
undoubtedly a much larger thermal boundary layer thickness than could be 
accommodated by the Hiemenz flow, but quasi-periodic nonlinear removal of 
heated matter from the region appears as an additional mechanism that could 
help to account for the observed effects as conjectured by Hodson and 
Nagib (1975) — see Section 6a and Figure 8. 

If the dimensionless difference of temperature from its wall value is 

AT = {T - T(ri) }/ (T - T ) and its linear perturbation is 0, (n) > the 
CU CO * 

differential equation for the latter, in the Hiemenz flow, reads: 

e^+ Pr$Q k^0^ = -PrAT* 

where Pr is the Prandtl number. Comparing this with Equation (13) for 
, we readily identify the thermal diffusion in n and C directions 
x 2 

as 6^ and -k 6^^ and the convective derivative -Vq 3G/9ri as Pr'I’Q 6| . 

The source-like term on the right represents modifications of the otherwise 
stratified temperature by the cellular "mixing motion" and plays the 

role of the stretching term -(9 Vq/8ti)o 3^ = • o)^ in Equation (13). Both 

-AT'V-i and its aforementioned counterpart in the U-i equation, 

, vanish outside of the boundary layer whereas the stretching 
term continues to infinity as in the Hiemenz flow. So like , 

is decoupled from the inviscid amplified vortical motion outside of the 

boundary layer and for Prandtl numbers near unity, i.e. for most gases, has 
n profiles very similar to that of in Figure 17. Typical profiles of 

linearized temperature perturbations can be found in Sutera (1965), Figure 9, 
and Williams (1968), Figure 16. The isotherms of the combined mean and 
perturbation field are essentially cosine like with higher temperatures 
occurring at z = 0 and X , where the perturbed flow is away from the 
wall (i.e. something like the V curve itself in lower Figure 20, but less 
exaggerated) . The isotherms are straight at the wall because of the constancy 
of temperature imposed by the authors as boundary condition and again past T) of 
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4, the edge of the thermal diffusion from the wall in the V counter flow. 

They are the waviest near 71-1.2 — 1.4 , the maximum of the 0^ profile. 

Computer Differences? Guideline No, 4 

In the preceding section we relied heavily on computers to elicit the 
general properties and features of the velocity, pressure, and temperature 
fields.. But computing is still an art rather than science in complex 
Instability problems. Could differences in computational details and errors 
due to approximations be also amplified? Comparison of the results of 
Sutera (1965), his Figure 4, and Inger (1974), our Figure 17, indicate 
substantial differences in relative variations for the one case computed by 
Independent programs: k = 2/3 . (Inger actually used a compressible code, 

but the results in Figure 17 are for AT = 0 , Ap = Ay = 0 , and yield the 
same Hiemenz base flow as a check.) Sutera’ s maximum of and zero of 

(vortex center height) are located at r) - 2.4 , while Inger ’s are at 
rj 3.0 . From its maximum at ri=2.4 to ri = 8.5, Sutera 's drops 

by a factor of 5 while Inger ’s drops only by a factor of 2 between 

r| = 3 and 9 . The ratio of the negative inner maximum of to its 

positive outer maximum (Figure 17 ) is approximately 2.7 for Sutera and 4.6+ for 
Inger. It seems that we can trust the general features discussed in the 
preceding section but not necessarily the magnitudes. 


It is not clear whether the differences creep in from error amplification, 
from handling of the conditions at ’’computer infinity" or from other causes. 


The Condition at Infinity: Reprise 


To start his integration toward the wall from a finite r) , Inger uses 
as 3 nnptotic solutions which insure Vj^(°°) = 0 . Sutera like all the Brown 


University investigators formally asks that 


Vj^(-) 


V” (00) etc. 


0 as n 


and speaks only of boundedness of at large r) not of vanishing of . 

The details of his procedure are not clear nor is it clear v^hether he agrees 
with the analysis of Hammerlin (1955) which indicates that for k< 1 , 

V^(co) -i- 0 implies Vj^(«>) -*■ 0 as Gortler (1955) pointed out — see Section 5h. 


All these investigations deal with Hiemenz base flow. For non-Hleraenz outer 

flow, V' («») 0 would not necessarily imply V-. (®°) -^0 as mentioned in 

J ^ 

Section 5. 


Tani and lida, as reported by Tani (1974), make an explicit statement that 
the condition 0 as r|-»- is incompatible with their differential equations 

which generalize the Gortler-Hammerlin system in several respects. First a 
circular cylinder is considered (with allowance for the wake) by introducing 

cylindrical coordinates and obtaining the factor (l+/v/a q/R)”^ , which 
modifies the limit of Vq as U , see Equation (5) in Section 2 d. Since 

for large n , 0.648 , “ 1.765 Voo (rather than «» or ) 

as Tani and lida also account for the second-order boundary-layer 
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effects of nonzero curvature and for displacement thickness. The detailed 
equations are not shown in Tani (1974) but it is likely that the first effect 
Is the governing one as regards the approach to infinity. The same modifying 
coordinate factor of Equation (9) is present in the definitions of the 
perturbations and . 

Flow Fields for Decreasing k Hore Computer Differences? 

Tani and lida's 1974 non*^Liemenzian system included a mixtiare of first- 
order and higher-order terms, which presumably caused (°o) to be a nonvanish- 
ing constant dependent on the wave number The smallest occurred for 

k=0, 3 5, corresponding to a longer wavelength! = 2 tt/ 0,35 A>/a = 27.70 6* ~ 7.49(5. 
Thus each vortex cell is just less than four Hiemenz boundary layer thicknesses 
wide. The corresponding n profiles in Figure 18 indicate the location of 
the linearized vortex center to be near n of 4.5 or almost two boundary- 
layer thicknesses out. It is curious that Inger’s adiabatic solutions with 
Hiemenz base flow (which have V 2 (<») 0 ) for the sequence k = 1.0, 0.667 , 

0.5, 0.2 show no outward shift of the vortex centers. The one striking 

trend in Inger’s sequence is the strong relative decrease of with 

respect to as k decreases. In fact his vanishes for all practical 

purposes at four Hiemenz boundary layer thicknesses from the wall for 
k *= 0.2: according to him, longer wavelengths flatten out the motion and 

inhibit the normal velocity perturbations. Tani and lida’s solutions in 
Figure 18 ( V]^(<») = min ^ 0 , curvature taken into account, etc.) appear then 

quite different, in particular in regard to the vigorous motion relative 

to the motion for such a relatively long wavelength. 

lida (1978a) systematically sorted out the orders of magnitudes in the 
rpani— lida model in preparation for his weakly nonlinear treatment (1978b) . To 
the first order (O(A-C^)^ in lida's notation)^ Hammerlin results of Section 5h 
were recovered. The only first-order profile printed in lida's paper is for 
k = 0,5, so that a direct comparison with Fig 18 to observe the differences is 
not possible. However, comparison with Inger's 1974 k = 0,5 profiles reveals 
that lida's results do not follow Inger's trend described above, In fact, 
lida's ratio V^max/W^max exceeds two, just as in Fig 18, whereas Inger's falls 
below 0,5^ 


The Hainzl-Hammerlin Contradiction and Possible Resolution 

One more clue for the far— field near— field relationship, ^ = 0 , should be 
mentioned, especially since it seems to have been forgotten. Gortler was 
well aware of the basic difficulties brought on by the parameter less , unbounded 
Hiemenz base flow. At his instigation Hainzl (1965) turned to a circular cylin- 
der and first found "a new rigorous solution of the steady Navier-Stokes equa- 
tions as a replacement for the Hiemenz flow”. Physically, Hainzl 's new base flow, 
see Fig 19, represents an incompressible field induc ed b y applying infinite suc- 
tion along the rear stagnation line . If the n = y magnification were 

applied, this flow could not be distinguished from Hiemenz', in rectangular 
coordinates fitted to the front stagnation line, and presumably from that of Tani 
and lida in cylindrical coordinates for which is given by Equation (9) . In 
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1968, Gortler privately described the situation as puzzling. 


One can surmise that the puzzle centers on the asymptotic behavior of the 
stream function in the limits r)~*oo and y oo . Neglecting the displace-^ 
ment effect for the limitin g pro cesses to follow Hainzl's stream function 
(Equation (5) ) for y> 4 /^/ a ^ i.e. beyond the boundary layer thickness S ,is 

#^(►7) = -in(l +-|) = /n(l + . (IS) 

For y values large with respect to 5 but small with respect to R , 

y/R = /v/a T\/R ~ of Hiemenz. The corresponding expression for the 
base velocity Vq beyond the boundary layer is according to Equation (9) 

Vqjj = -*^av jiln(l+ y/R)j(l+ y/R) ^ which in the intermediate limit 
^ ^*q(ti)( 1+ A7/a n/R) ^ ^ the Tani-Iida expression. As ri grows large, 

the latter approaches -aR , not , when the Hiemenz approximation is 
Imbedded in a potential flow with uniform upstream V^o > tiot zero, the 
correct y limit of Vq^ when the Hiemenz flow is imbedded in Hainzl’s new 

field solution. In the light of our discussion in Section 5c of the vorticity 
amplification source BVq/Sy there is ample opportunity for differences 

in results depending on which limit is taken and how carefully any matching 
of functions is performed. In Hainzl's case, vorticity disturbances oncoming 
from large y are first inviscidly depressed (cf narrowing of stream tubes 
coming from far upstream in Figure 19) as well as depleted by diffusion. Only 
for y/R of less than 0.5 or so (in the 0 wedge of convergence of the 
perturbation equations) is there inviscid amplification. How far does 17 -»• oo 

really extend when the Hiemenz flow is analytically contained in a given 
external flow, such as that of Hainzl, and one is using a computer? The 
puzzle could be explained if Hammerlin (1955) worked within the limit of the 
inviscidly amplifying region near the cylinder and Hainzl apparently beyond 
it for his base flow. 
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Figure 1 . 


Hiemenz two-dimensional stagnation flow. (After Schlichting,1968 . ) 
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Figure 2, Examples of imbedding of Hiemenz flow in outer irrotational flows 
which provide a characteristic length. 
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SIDE VIEW 


Figure 3 . 


Schematic of counterflow vortex formation due to impingement of 
steady wakes of "rods" on stagnation region of bluff bodies. 
(Courtesy of Hodson and Nagib# 1975.) 






Figure 4. Centerplane smoke visualization of three outer streamlines in counterflow vortex forma- 
tions caused by placing a prismatic obstacle in a laminar wall boundary layer, (Cour- 
vo tesy of R. Norman, 1972.) 



d=.69mm D= 23.75mm 



(b) Re^=365 ReQ = 2730 Ay/d = 32 


d= 3.18 mm D = 23.75 mm 


Figure 5. Side-view dye visualization of vortex formation due to a steady and 
a mildly unsteady single rod wakes . In (b) far— upstream laminar 
dye injection causes small momentum defects at the edges of the wake which may 
influence near-singular flow conditions. (Courtesy of Hodson and Nagib, 1975.) 
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Figure 6, 


Empirical correlation of threshold conditions for counterflow vor- 
tex formation for single and multiple wakes impinging on a rectang 
ular cylinder, (Courtesy of Nagib, Hodson, and Roadman, 1977.) 
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side-view dye visualization of strong unsteady interaction between 
a counterflow vortex pair and a "natural" perturbation P of the 
vorticity in the rod wake at Re^ = 90, (From a movie of Hodson 
and Nagib, 1975.) 
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Figure 9. Oscillographs of u'(t), the perturbation velocity perpendicular 
to a hot wire located at various normal distances y from the 
cylinder, presumably in the stagnation plane. Water tank experiments: 

Voo = 15 cm/s; S = 1.8 mm; time = 0,02 sec/div; sensitivity lOOmV/div; 

Tu /«s/ ,85% generated by a screen riding 20 cm upstream of the cylinder; Tu 
measured 10 cm downstream^ and 10 cm above the stagnation plane; screen mesh: 
12/cm; wire diameter: 0.31 mm, (Courtesy of H. Hassler, 1971, and private 
communication , ) 
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Figure 10. Oscillographs of u' (t) , the perturbation velocity perpendicular 
to a hot wire located at a normal distance y = 2 mm from the 
cylinder in planes at angles 'f from the stagnation plane . Conditions as in 
Figure 9. (Courtesy of H. Hassler, 1971.) 


97 


I 


00 



Figure 11, Vertical view downward on the 12 cm diameter cylinder moving in a water tank with 

Voo ^ 11.9 on/s into screen-generated turbulence, Tu, of approximately 0,6%. Hydrogen- 
bubble time lines are formed at a horizontal cathode-pulsed wire location at y = 2mm, = 11° at the 
top of the picture. The time lines are progressively deformed by the upwelling motion of cin irregular 
series of imsteady vortex pairs. Originally, the bubbles move towards the plane of symmetry of each 
pair, are uplifted and then thin out as they move away from the plane of symmetry. Ultimately the 
bubbles tend to concentrate at the core of the vortices. The time-line shapes near "f = 80 correspond 
to the integrated history of the motion, not to the local vortex strengths. (Courtesy of P, Colak- 
Antic, 1971, and private communication.) 




Figure 12. Strengthening of vortices with increase in V.* to 19.4 cm/s, other conditions remaining 

. . ' . as in Figure 11. Stereoscopic viewing identifies the seemingly closed spirals at left 

•:'y center as a succession of time-line tongues stretched by rotation into mushroom shapes 

with low bubble density at the top. The concentration of bubbles in the cores develops 
somewhat irregularly. Stronger vortex pairs often have weak neighbors and there is some amalgamation, 
(Courtesy of P, 6olak-Antic, 1971, and private communication.) 
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Figure 13. Distribution of u(y,z)^ the short-tjLme averages (over one second) 
of velocity perpendicular to a hot wire located at spanwise loca- 
tions z and normal distances y in the plane = 4 . Other conditions as in 
Figure 9. (Courtesy of H. Hassler, 1971.) 
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Figure 14. Distribution of one-second averages of velocity u' perpendicular 
to a hot wire at normal distances y from the cylinder, at various 
•f angles and a fixed z location. Other conditions as in Figure 9. (Courtesy 
of H. Hassler, 1971.) 
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Figure 15, Distribution of one-second averages of velocity u(2,Zf^) perpen- 
dicular to a hot wire located at spanwise locations z, a fixed 
normal distance y = 2 mm, in planes at angles from the stagnation plane. 
Other conditions as in Figure 9. (Courtesy of H. Hassler, 1971.) 
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Figure 16. Distribution of one-second averages of velocity u(y,z) perpen- 
dicular to a hot wire located at spanwise locations z and normal 
distances y in the stagnation plane 0 , Other conditions as in Figure 9, 

(Courtesy of H. Hassler, 1971.) 
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Figure 17, Steady disturbance field profiles of velocity components, static 
pressure and vorticity; eigenvalues for prescribed wave number 
k = 2/3, (Courtesy of G, Inger, 1974.) 




Figiare 18, 


Velocity and vorticity solutions of the steady (neutral) 

^ eigenvalue problem, subject to asymptotic condition that 

be a mxnimum. This condition selected k = 0.35. (After Tani, 1974.) 
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FULL NON-LINEAR STREAMFIELD 



CORRESPONDING LINEAR PERTURBATIONS 

Figure 20, Second-order streamline pattern (Courtesy of Kestin and Wood, 
1970) and corresponding linear disturbance fields. 
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Figure 21. Nondimens ionali zed, experimentally observed, time-averaged span- 

wise wavelengths in flows with various turbulence levels as 
functions of cylinder Reynolds number, (After Tani , 1974.) Diamond-shaped 
data is from Figiire 16 and is probably unreliable because of probe interfer- 
ence . 


106 




2 . Govtmmtnt Accanion No. 


3. Racipiant'i Catalog No. 


1. Report No. 

NASA CR-3231 


4. Title and Subtitle 

S. Report Date 

ON THE QUESTION OF INSTABILITIES 

December 1979 

UPSTREAM OF CYLINDRICAL BODIES 

6. Performing Organization Code 

7. Author(s) 

8. Performing Orgenization Report No. 


I.I.T. Fluids and Heat 

Mark V. Morkovin 

Transfer Report R79-3 


10. Work Unit No. 

9. Performing Organization Name and Addreu 

Illinois Institute of Technology 
MMAE Department 


11. Contract or Grant No. 

Chicago, IL. 60616 

NSG-1120 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Contractor Report 

National Aeronautics and Space Administration 

11/15/74 - 4/30/77 

Washington, D.C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


Langley Technical Monitor: Ivan E. Beckwith 


Final report 



16. Abstract 


In an attempt to understand the unsteady vortical phenomena in perturbed stagnation 
regions of cylindrical bodies , a critical review of the theoretical and experimental 
evidence was made. Current theory is revealed to be incomplete, incorrect, or 
inapplicable to the phenomena observed experimentally- The formalistic approach 
via the "principle" of exchange of instabilities should most likely be replaced 
by a "forced-disturbance" approach. Also, many false conclusions were reached by 
ignoring that treatment of the base and perturbed flows in Hiemenz coordinate n 
is asymptotic in nature. Almost surely the techniques of matched asymptotic 
expansions will have to be used to capture correctly the diffusive and vorticity 
amplifying processes of the disturbances vis-a-vis the mean-flow boundary layer 
and outer potential field as ri and y/diameter approach infinity. There are also 
serious uncertainties in the experiments which are discussed in detail. 






17. Key Words (Suggested by Author(s)) 

18. Distribution Statement 


Instability Theory; Stagnation flow; 

Flow about a bluff body; Vorticity 
amplification; Turbulent boundary layer; 
Laminar boundary layer; Enhanced heat 
transfer. 


Unclassified - Unlimited 


Subject Category 34 


19. Security Oaaif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

110 

$6.50 


For sale by the National Technical Information Service. Springfield, Virginia 22161 


NASA-Langley. 1979 





















